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Method and system for estimating the position of a mobile device. 



FIELD OF THE INVENTION 

5 

The present invention relates to a location metliod and in particular but not 
exclusively to a method of determining the location of a mobile station In a wireless 
communications network. 

10 BACKGROUND OF THE INVENTION 

Wireless cellular communication networks are known. In these networks, the area 
covered by the network is divided Into a number of cells. Each of the cells has 
associated with it a base transceiver station. The base transceiver stations are 
15 arranged to communicate with mobile devices located in the cells. The mobile 
devices can take any suitable form and are typically mobile telephones. 

The need for efficient accurate positioning of the mobile telephones has increased. 
In the USA, network operators must be able to provide the location of mobile 

20 subscribers making emergency calls. Similar proposals are currently being 
considered in Europe. Additionally, commercial services such as tracking services 
(that is the emergency service mentioned above, people locating, valuable assets 
location etc), finding/guidance services (proximity services such as yellow pages, 
direction indicators, point of interest locators etc) and notification services (targeted 

25 advertisements, traffic alerts, weather alerts, bus/train location, guided tours etc) are 
currently being proposed. 

In the GSM (global system for mobile communications) standard, four positioning 
methods have been included: cell identity and timing advance, time of arrival, 
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enhanced observed time difference (E-OTD) and a method based on GPS (global 
positioning system technology either as standalone GPS or assisted GPS). 

The time of arrival method is able to locate handsets with standard software but 
5 requires the installation of new network elements such as location measurement 
units at every base station site. The enhanced - observed time difference method 
requires the installation of location management units at every two to five base 
transceiver stations and a software modification in the handset. The assisted GPS 
method requires installation of a GPS receiver and possibly also location 
10 measurement units besides the integration of a GPS receiver into a handset. All of 
these methods require the introduction of a new network element or corresponding 
functionality responsible for location calculation called the serving mobile location 
centre SMLC. 

15 The timely deployment of location services for all users including those having 
handsets which do not include the necessary software or hardware, requires that 
measurements already available in cellular networks should be used. These 
techniques are important as they allow operators and service providers to start 
offering location based services to all customers with minimal additional costs whilst 

20 waiting for more accurate and sophisticated location technologies to be available. 
From the technical point of view, even when technology such as enhanced-observed 
time difference and assisted GPS are fully available, network based software 
solutions will still be needed as backup methods when the new standardised 
solutions fall or when the requested accuracy can be met with such a method. 

25 Network based software techniques can also be used as an initial guess for the 
algorithms used to Implement one of the standard solutions in order to improve the 
accuracy or speed of convergence of those algorithms. 



30 



SUMMARY OF THE INVENTION 
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According to a first aspect of tlie present invention, there is provided a metliod of 
estimating the location of a mobile device, comprising the steps of: collecting location 
Infomnation; selecting at least one of a plurality of different location methods to 
provide a location estimate; and providing a location estimate based on the at least 
5 one selected location method. 

According to a second aspect of the present Invention, there is provided a system for 
estimating the location of a mobile device, comprising: means for collecting location 
information; means for selecting at least one of a plurality of different location 
10 methods to provide a location estimate; and means for providing a location estimate 
based on the at least one selected location method. 

BRIEF DESCRIPTION OF DRAWINGS 

15 For a better understanding of the present invention and as to how the same may be 
carried into effect, reference will now be made by way of example only to the 
accompanying drawings in which: 

Figure 1 shows wireless cellular networl< to which embodiments of the present 
Invention can be applied. 
20 Figure 2 shows a mobile device served by three base transceiver stations in a 
wireless cellular network as shown in figure 1 . 

Figure 3 shows a geometric representation of a base transceiver station and mobile 
device In a wireless cellular network as shown in figure 1 . 
Figure 4 shows an example of antenna gain in a transceiver as shown in figure 1 . 
25 Figure 5 shows another example of antenna gain in a transceiver as shown in figure 
4. 

Figure 6 shows approximated radiation patterns as found in transceivers as shown in 
figure 1. 
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Figure 7 shows the Okumura-Hata path loss graphs for modelled wireless cellular 
networks as shown in figure 1 . 

Figure 8 shows a geometric estimate for the confidence region of a location estimate 
for a serving cell In a wireless network as shown in figure 1. 
5 Figure 9 shows the confidence region of a location estimate provided by a series of 
cell's cell identity. 

Figure 10 shows a geometric representation of the coverage of a cell In a wireless 
cellular network as shown in figure 1 . 

Figure 1 1 shows a series of plots of the probability density function of a location 
10 estimate for various values of path loss/attenuation In a wireless cellular network as 
shown in figure 1 . 

Figure 12 shows a plot of ho against path loss/attenuation In a wireless cellular 
network as shown in figure 1 . 

Figure 13 shows an example of the serving area of a cell In a wireless cellular 
15 network as shown in figure 1 . 

Figure 14 shows the Circular crown confidence region for a location estimate 
produced from a constant TA value. 

Figure 15 shows a further geometric view of the confidence region for a location 
estimate as shown In figure 14. 
20 Figure 16 shows the geometry used to calculate the confidence region in a CI 
location estimate. 

Figure 17 shows a flow diagram detailing the steps to provide location estimates in 
embodiments of the invention. 

25 

DETAILED DESCRIPTION OF EMBODIMENTS OF THE PRESENT INVENTION 

Reference is made to Figure 1 which shows schematically a wireless 
telecommunications network 2 to which embodiments of the present invention can be 
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applied. The area covered by the network is divided up Into cells 4, Associated with 
each cell is a base station 6. In other embodiments of the invention base stations 6 
may be typically referred to as base transceiver stations (BTS). The base stations 6 
are arranged to communicate with user equipment 8 via a wireless connection. The 
5 user equipment is typically a mobile device such as a mobile telephone, a computer, 
a personal digital assistant (PDA) or the like. In other embodiments of the invention 
the user equipment may also be typically referred to as a mobile station (MS). 

Embodiments of the present invention are arranged to estimate the location of a 
10 mobile device (IVIS) and an associated confidence region. The confidence region is 
the region in which the mobile device (IVIS) can be expected to be located with a 
given probability. In other words, an estimated position of the mobile device (MS) can 
be calculated but that estimate will not be 100% accurate. The confidence region 
defines an area in which it is possible to be reasonably certain that the mobile device 
15 (MS) is located. 

Embodiments of the present invention are arranged to combine a number of location 
algorithms which will be described in more detail hereinafter. 

20 The location algorithms used in embodiments of the present invention use cell 
identity (CI), timing advance (TA) and received signal strength (RX) measurements. 
Preferred embodiments use three types of algorithms: 

1) algorithms based on the cell identity. These are referred to as CI algorithms; 

2) algorithms based on the cell identity and timing advance information. These are 
25 referred to as CI +TA algorithms; and 

3) algorithms based on the cell identity and received signal strength measurements. 
These are referred to as CI + RX algorithms. 
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These different algorithms will be described in more detail hereinafter. For 
convenience these algorithms are divided into two categories. The first includes the 
CI + RX algorithms and the second includes the other two types of algorithm, that is 
CI and CI + TA algorithms. 

5 

It should be appreciated that in some networks, for example the GSM network, the 
cell identity, timing advance and received signal strength data is already available 
and used for other purposes. This means that at least some of the algorithms can be 
implemented without requiring any modification to the existing handsets. 

10 

The location algorithms used in embodiments of the present invention require a 
certain set of simultaneous equations to be solved. The simultaneous equations, 
depending on the statistical assumptions of the measurements can be non-linear or 
linear in the unknowns. If the equations are non-linear an iterative method has to be 
15 adopted to find a solution while if the simultaneous equations are linear, a solution in 
closed form exists. Closed form algorithms are computationally lighter than iterative 
algorithms. 

CI + RX algorithms 

20 

The class of algorithms based on the cell identity and received signal strength will 
now be described. 

Reference is made to Figure 2 which shows the principles of location estimates 
25 based on received signal strength. Received signal levels are measurements of the 
levels of received signals from the serving base transceiver station (BTS) 6 (the base 
transceiver station (BTS) with which the mobile device (IVIS) 8 is associated) and the 
six strongest neighbours. Of course other numbers of base transceiver stations can 
be used. These measurements are performed by the mobile device (MS) 8 and 
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reported to the fixed part of the network when the mobile device (MS) 8 is in a 
dedicated mode, in the idle mode, the mobile device 10 measures the received 
signal level from the best server, that is the base transceiver station (BTS) on which it 
is camped and the six strongest neighbours. However, these measurements cannot 
be reported to the fixed part of the network as there is no connection with any of the 
base transceiver stations. 

The level received by a mobile device or more precisely attenuation that the received 
signal has experienced, depends on the reciprocal position of the mobile device and 
base transceiver station involved. The received signal level from multiple base 
transceiver stations can then be combined to estimate the location of the mobile 
device. 

Embodiments of the present invention belong to the class of network based software 
solutions (NBSS). In particular, embodiments of the present invention described in 
this section are based on the use of signal level measurements and focus on 
applications where coverage prediction maps are not available. However, 
embodiments of the present invention can be used where coverage prediction maps 
are available. Coverage prediction maps are used In some methods where the 
mobile device location is estimated as the location on a coverage map where the 
values of predicted signal strengths best match the levels actually measured. 

The measurement/lnfomnation needed to implement the algorithms described 
hereinafter can be separated into networic configuration parameters such as: base 
station coordinates; sector orientation and base transceiver antenna radiation 
pattems In the case of sectorised cells; and maximum base transceiver station 
downlink transmission power, and measurements such as the level of signals 
received by the mobile device. In addition to that data, knowledge of suitable models 
to link level measurements to the reciprocal position of the mobile devices to be 
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located and the base transceiver stations used in the positioning procedure are 
required. These will be discussed in more detail below. 

The location methods outlined hereinafter estimate the mobile device coordinates at 
the location where a certain function has its minimum value. The function Is obtained 
by combining level observations which are estimates of the attenuation experienced 
by the signals received by the mobile device and transmitted by the base transceiver 
stations Involved. Level observations are estimated by subtracting from the levels 
received by the mobile device, the contribution of the base transceiver stations 
antenna radiation patterns (which is dependent on the reciprocal angular position 
between the mobile device and the base transceiver stations), the path loss (which is 
dependent on the distance between the mobile device and the respective base 
transceiver station) and other constant factors such as the base transceiver stations 
transmission power, cable losses, antenna losses etc. 

Clearly, the accuracy of the results provided by the algorithm depends on the 
accuracy of the information available. The accuracy of the result can be improved by 
including accurate definitions of the antenna radiation pattems and finely tuned 
propagation models. These models. In principle, do not need to be expressed in an 
analytical fomi but may be included In the location algorithms as look-up tables from 
which, given a certain mobile device to base transceiver station angle of arrival, the 
antenna gains can be retrieved. However, this further refinement can be omitted 
from embodiments of the present invention. 

It is preferred that the networl< information such as base transceiver station 
coordinates, sector infomnation, antenna radiation patterns and other adjustable 
parameters should be constantly updated. 
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In preferred embodiments of the present invention, the models used to represent 
signal level measurements may be continuously adaptive to changes in the radio 
environment, preferably in an automatic way. One way to accomplish this is to adjust 
the models based on statistical comparisons between a certain set of collected 
measurements and the corresponding true quantities. This however requires the 
exact position of the mobile device to be known in order to calculate the true 
quantities and to compare them with the measurements. When the measurement 
models are derived offline, drive tests can be performed collecting at the same time 
level measurements and mobile device location information by means of a GPS 
(Global Positioning System) receiver co-located with the mobile device. There may 
be other techniques which are implemented which provide a more accurate location 
technique, for example, those using an estimated observe time difference or arrival 
technique or an assisted GPS location estimate, that allow the determination of the 
measurement models on-line. 

The algorithms which are described hereinafter are applications of the maximum 
likelihood principle to estimate the mobile device coordinates by processing a set of 
observations from the level measurements. Different maximum likelihood 
approaches are set out. 

The location algorithms use signal level measurements to estimate the mobile device 
coordinates, x and y, along with the variants of the slow fading assumed to be 
equal for all the base stations involved. 

As explained above In embodiments of the present Invention described in this section 
the input data required by the algorithms includes signal strength measurements 
collected by the mobile device (MS), network parameters, and a suitable path-loss 
law - in other words a suitable model for the wireless transmission characteristics. 
The list of Input data follows. 
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5 



Average power of the signals received by the MS from N base transceiver stations 
(BTS's) measured In decibels, PV: 

(1) 

Coordinates of the N measured BTS's in meters, (x'.y"): 

(^W) (^ = 1 N) ;H = M=m 

Orientation of the sectors of the N base transceiver stations (BTS's) (in the case of 
10 sectorized cells) in radians measured counter-clockwise from the x direction, <|>'b: 

Maximum radiated power from N measured BTS's In decibels: 

=Pt+ Gi,moar + Gr.max - LoSSeS* J [i^,^] = dB 

15 

The maximum radiated power, P't.max. represents the maximum power at the output 
of the i-th BTS antenna in the direction of maximum gain. It includes the transmitted 
power, P'l, the maximum gain of the BTS transmit antenna, G't.max. antennas losses, 
cable losses, etc.; all measured in dB. To simplify the mathematical model, P't.niax 
20 includes also the maximum gain of the MS receive antenna, Gr.max- 

Combined Transmit-Receive antenna pattern for the i-th BTS (i = 1 N) in decibels: 

APl(^'{x,y))^APi{^l;'{x,y)-<f>^)+APr{^'{x,^ ; [AF^]=dB 

25 The combined radiation pattern, AP'tr(V(x,y)), represents the gain Introduced by the 
antenna Installed at the BTS site and the antenna Installed In the handset. The gain 
depends on the reciprocal orientation of such antennas. In the embodiment 
represented by the above equation, <|>m Is the orientation of the antenna Installed at 
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25 



the MS, (|>'b is the orientation of the antenna installed at the i-th BTS and V(x,y) is the 
angle of arrival of the signal transmitted by the i-th BTS and received by the mobile 
device (MS), measured in radians counterclockwise fix>m the x direction (see figure 
3): 



i>'ix,y) = tan-^ ^1-1 ; U*! = rad 
X' — a; "-J 



(6) 



The combined radiation pattern APV, is a function of the MS coordinates (x,y). It 
includes the radiation pattern of the transmit antenna installed at the i-th BTS, 
AP't(e), and the radiation pattern of the antenna Installed at the mobile device 
10 (MS), APr(e), In embodiments of the invention the antenna at the MS is omni- 
directional and the orientation of such antenna is unknown; thus APr(6)= OdB and 
<j>M = 0 rad. 

The path-loss law for the propagation between MS and the i-th BTS measured in 
15 decibels: 



PL»((f(a;,y)) ; [PL»] = dB 



The path-loss PL'(d'(x,y)) represents the attenuation experienced by the signal 
transmitted by the i-th BTS as it propagates further away from the transmit 
20 antenna. It is expressed as a function of the distance between the MS and the i- 
th BTS, d', which in turn depends on the MS coordinates: 

d^M = yfi^i - a;)2 + (y^ - 



(8) 



The output of the location algorithms comprises: 
An Estimate of the MS location in meters: 

(9) 
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An Estimate of the variance of the log-normal slow fading in decibeis: 

(10) 

5 As mentioned previously, the algorithms apply maximum likelihood estimation 
principles. They estimate the mobile device (MS) location and variation of the slow 
fading by minimising certain scalar functions. The functions are obtained by 
combining level observations which are measurements of the attenuation 
experienced by the signal received by the mobile device (MS) with the corresponding 

10 expected quantities, determined as a combination of the base transceiver stations 
antenna radiation patterns (which are dependent on the reciprocal angular position 
between the mobile station and respective base transceiver stations) and path loss 
(which is dependent on the distance between the mobile station and the base 
transceiver stations). By minimising a certain cost function, (a cost function is a 

15 generally applied term in mathematics for the description of optimization problems), 
the algorithms find the value of the unknown parameters which globally minimise the 
difference between the observed attenuation and the expected attenuation. The cost 
function is a measure of the expected enters in the estimation of the location of the 
mobile device (MS). Therefore by minimising the cost function the expected location 

20 error is also minimised. Embodiments of the invention use the following algorithms 
depending upon assumed characteristics of the wireless environment. 

Embodiments of the present invention perform the first algorithm A where the 
wireless environment is assumed to comprise correlated slow fading 
25 characteristics and where the slow fading is assumed to have equal variance 
statistics. In other words the signal transmitted by each base transceiver station 
(BTS) to the mobile device (MS) has a similar but not identical characteristics. 
This type of assumption is accurate where the transmission paths between base 
transceiver stations (BTS) and mobile device are similar. 



wo 2004/023155 



PCT/IB2002/003644 



Algorithm A: Maximum Likelihood Estimation with Level Observations Assuming 
Correlated Slow Fading and Equal Variance 

5 1. Calculate the l-th level observation, U, by subtracting from the l-th 
measured received power, pV, the maximum power radiated by the i-th BTS, 



10 L* is the total attenuation experienced by the signal transmitted by the i-th BTS 
while propagating toward the IVIS. The total attenuation depends on path-loss, 
gains introduced by BTS antenna and MS antenna, fluctuations of the radio 
channel, etc. 

15 2. Stack the level observations from N BTS's in vector L: 



r« = p« _ pi 



(11) 





(12) 



3. Solve the minimization problem: 




L y 



(13) 



20 



where the cost function F (x,y; <7u^) Is defined as follows: 

F (x,y; or„2) = Ina^^ + In \Ti,(x,y)\ + inj,{x,y)f vi,-^{x,y) [L - mL(x,y)] 



(14) 



and 
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AM = -PL' - API {^\x,y)) 



(15) 



(16) 



M = 1,...,JV 



(17) 



Pu'''(x,y) is the cross-correlation of the slow fading affecting the signals 
propagating from BTS' and BTS^ toward the MS. 

Embodiments of the present invention perform the second algorithm B where the 
wireless environment is assumed to comprise uncorrelated slow fading 
characteristics and where the slow fading is assumed to have equal variance 
statistics. In other words the signal transmitted by each base transceiver station 
(BTS) to the mobile device (MS) has unrelated characteristics. This type of 
assumption is accurate where the transmission paths between base transceiver 
stations (BTS) and mobile device (MS) have no similar components. 

Algorithm 8: Maximum Likelihood Estimation with Level Observations Assuming 
Uncorrelated Slow Fading and Equal Variance 

1. Calculate the i-th level observation by subtracting from the l-th measured 
received power, P'r, the maximum power radiated by the i-th BTS, P't,inax: 



2. Stack level observations from N BTS's in vector L: 




(19) 



3. Solve the minimization problem: 
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it] 



= axg imn F (a;,y) 




(20) 



where the cost function F(x,y) is defined as follows: 



N 




2 



(21) 



and Dxy Is the domain of existence of x and y. Several possible definitions for Dxy are 
given later in tiie present application. 

4. Calculate al as 



Algorithm B differs from Algorithm A In the definition of the cost function. At the basis 
of this different definition Is a different model for the slow fading affecting the signals 
transmitted by two BTS's. In Algorithm A the fading is assumed correlated while In 
Algorithm B is assumed uncorrelated; this results in a simpler definition of the cost 
function of Algorithm B. 

Embodiments of the present invention perform the third algorithm C where the 
wireless environment is assumed to comprise uncorrelated slow fading 
characteristics and where the slow fading is assumed to have equal variance 
statistics, in other words the signal transmitted by each base transceiver station 
(BTS) to the mobile device (MS) has unrelated characteristics. This type of 
assumption is accurate where the transmission paths between base transceiver 
stations (BTS) and mobile device (IVIS) have no similar components. 



(22) 
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Algorithm C: Maximum Likelihood Estimation with Level Difference Observations 
Assuming Uncon-eiated Slow Fading and Equal Variance 

1. Calculate the i-th level observation by subtracting from the l-th measured 
5 received power, Pt, the maximum power radiated by the l-th BTS, P't,, 

L' = I^-Pl^ I i = l,...,N 



:,max- 



(23) 

2. Calculate the j-th level difference observation by subtracting the j-th level 
observation from the level observation L^ taken as reference: 



10 Di = L^-U ■ j^2,...,N 



(24) 



3. Stack the N - 1 difference of level observations in a vector D: 

°=[^' (25, 

15 4. Solve the minimization problem 



[ ^ ] = axg imn F {x,y) 

y J 



(26) 

where 

i=2 ■'V \j=2 J (27) 

20 

and 

l^(^,y) = - [PL^ {d\x,y)) - PU {(Pix,y))] - [API {i,\x,y)) - API (V^(x,y))] 



(28) 
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Dxy is the domain of existence of x and y. . Several possible definitions for selecting 
Dxy are given later, 

5. Calculate 6l as 

5 (fu^ = F (x,y) (29) 

Algorithm C differs from Algorithms A and B in the definition of the observations 
and, as a consequence, in the definition of the cost function. The observations 
considered in Algorithm C are differences in the attenuation experienced by 
10 signals transmitted by two different BTS's and received by the MS. To slmpllly 
Algorithm C, moreover, the slow fading processes affecting signals transmitted by 
two BTS's are assumed uncorrelated, analogously as in Algorithm B. 

Algorithms A, B and C described hereinbefore are non-linear algorithms. 
15 The above location algorithms are hereafter further detailed and described. 

Figure 3 shows the basic geometry of the problem. Figure 3 comprises a mobile 
device (MS) 8. and the l-th base transceiver station (BTS) 6. The MS 8 and BTS 6 
exist in a region defined by a Cartesian co-ordinate system 301. All angles are 

20 defined as being defined from the x-axIs in an anti-clockwise direction. The i-th 
BTS 6 is located at the point defined as (x'.y') and arranged to broadcast and 
receive with a maximum gain direction 303 on an angle defined as ^'b- The mobile 
device MS 8 is located at the point (x,y) and is arranged also to transmit and 
receive with a maximum gain direction 305 on an angle ^m- The MS 8 is 

25 positioned in relation to the i-th BTS 6 by a line 307 comprising length d'(x,y) and 
by an angle \|/'(x,y). In further embodiments of the invention alternative co-ordinate 
systems to the Cartesian system are used. In other embodiments of the invention 
the polar reference system centred on BTS' with the distance between MS and 
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BTS', d, defined as the radial coordinate and the Angle Of Amval of the signal 
received by the mobile device. \|/'. defined as the angular coordinate: 

The distance between the mobile device (IVIS) and the i-th base transceiver station 
5 (BTS) is defined according to equation (30): 

Angle of arrival of the signal transmitted by the I-th BTS and received by the 
mobile device (MS) is defined according to equation (31): 



10 



15 



25 



X* — X 

(31) 

The transformation between x,y (Cartesian) and cf, \|/' (radial) coordinates is 
completed by the following formulas: 



{ 



X = + cos ^' 
y = y* + ti'sinV'* 



(32) 



The location algorithms as described in the above embodiments estimate the MS 
location by processing certain level observations. Such observations are obtained 
from signal strength measurements perfomied by the MS. This section derives a 
model for level observations from a model of signal strength measurements 
20 performed by the MS. Such measurements are, by definition, estimates of the 
average power of the received signals. A general model to calculate the average 
power received by the mobile at a particular location, (x.y), from the i-th base station 
is the following (all quantities are in dB): 

K(x,y) = Pi + Gi{x,y) + C?j(x,2/) - PL*(a;,y) - Losses* H- u\x,y) dB 



(33) 
Where 
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P'r(x,y) Is the power of the signal received by the MS at location (x,y); 
P't Is the BTS transmitting powers 

G'r(x,y) Is the MS antenna gain In the direction of the i-th BTS; 

C3't(x,y) is the antenna gain of the I-th BTS In the direction of the MS; 

PL(x,y) is the path-loss determined by the propagation path between MS and BTS; 

the temi "Losses'" takes into account the losses due to antenna feeder, cables, 

duplex, divider, etc. 

u'(x,y) is the shadow fading affecting the signal transmitted by the I-th BTS. It Is 
generally modelled as a random variable with log-normal distribution (i.e, u'(x,y) 
measured In dB is a Gaussian random variable) with standard deviation a'u. Typical 
values for a'u range from 5 to 10 decibels: 

u*{^,y) - ^(0,4) [u'] = dB ^34j 

Equation (33) is the starting point to study techniques for locating mobiles by means 
of observations derived from observed level (or received level RXLEV) 
measurements. It is a general enough model to allow introduction of several 
parameters affecting the propagation such as antenna radiation patterns, path loss 
and random fluctuations. 

Antenna gains are usually expressed analytically as a function of an angle, 9, which 
describes the angular distribution of the power radiated by the antenna when it is 
connected to a transmitter. In case of directional antennas, the direction of maximum 
gain is identified by the direction 6 = 0: 

G{e = 0) = Grru^ ^35j 



The antenna gain G(e) can be separated in the sum of a constant term defining 
the maximum antenna gain, Gmax. and the radiation pattern, AP(e), which 
describes the angular distribution of the radiated power: 
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G{e) = Gmax - AP{e) ; AP{e)>OdB ; AP(5 = 0) = 0 dB 

(35) 

Figure 4 shows a graph of an antenna gain for an antenna as used in an 
embodiment of the invention. The graph of figure 4 comprises the y-axis 
representing antenna gain G(e) and the x-axis representing the orientation from 
the angle of maximum gain. The antenna gain Is symmetrical about the 9 = 0 
line. The antenna gain graph comprises a main lobe (or beam) 401 and four side 
lobes (or beams) 403, 405, 407, 409 either side of the main lobe. The main lobe 401 
comprises a maximum gain Gm centred at 6 = 0, The main lobe antenna gain 
decreases rapidly to zero either side of 9 = 0. The side lobes 403. 405, 407, 409 
comprise smaller maximum gains, and arranged so that the first side iobe 403 (the 
side lobe directly adjacent to the main lobe 401) has a larger maximum gain than the 
second side lobe 405 (the side lobe adjacent to the first side lobe 403 and the third 
side lobe 407). The second lobe 405, in turn, has a larger maximum gain than the 
third lobe 407 (the side lobe adjacent to the second side lobe 405 and the fourth side 
lobe 409), and the third side lobe 407 has a larger maximum gain than the fourth side 
lobe 409 (the side lobe adjacent to the third side lobe 407). 

When analytic expressions of the antennas- radiation patterns are not available 
an approximation is needed. 



The orientation of the BTS antenna, <|)'b (as shown in figure 3) specifies the 
direction toward which the maximum power Is radiated. This means that, given 
the (x.y) MS coordinates, M/'(x,y) (e.g., the angle In which BTS' "sees" the IVIS) 
can be calculated and the downlink transmission gain toward the IVIS can be 
expressed as follows: 
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In principle, if the MS Is equipped witli a directional antenna and the orientation 
of the antenna is known, an analogous formula can be used to describe the 
gain of the receive antenna at the MS side: 

5 

Although mobile devices can be equipped with directional antennas, it is not realistic to 
assume that the orientation of the MS antenna is known; thus in other embodiments 
of the present invention the MS antenna radiation pattern can be modelled with an 
averaged omni-directional pattern by imposing the following constraints 

The substitution of (36) and (37) in (33) results in the following compact 
expression for the received power: 

P; M = ^ API {^lj\x,y)) - ^V{x,y) + u\x,y) 



15 



(39) 

where the auxiliary definitions PVmax and AP'tr have been defined previously in 
equations (4) and (5) respectively. 



20 Several analytical expressions for the path-loss PU in equation (33) have been 
proposed in the scientific literature. In fact, this term is at the basis of the 
propagation loss prediction models based on which cellular operators design 
their networks. In order to calculate precisely the attenuation experienced by a 
signal travelling from a base station to a mobile, finely tuned prediction models 

25 including also information on topography and morphology of the environment 
should be considered. However, when terrain maps are not available, simplified 
models must be used. Such available propagation models are known in the art 
and in embodiments of the present invention the propagation model the path 
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loss as a function of the MS-to-BTS distance only with no dependence from the 
angle of arrival (AOA): 

With the above considerations in mind, the average power received by the 
mobile device (MS) at a particular location. (x,y), from the i-th base station pV 
can be expressed as follows: 



(41) 



where previous definitions for; the maximum radiated power from the i-th base 
transceiver station P't,^ax [equation (4)], the combined transmit-recelve antenna 
pattern APV [equation (5)], the path loss between the mobile device (MS) and 
the I-th base transceiver station (BTS) PL' [equation (7)] and the log-normal 
shadow fading affecting the signal transmitted by the I-th BTS u' [equation (34)] 
hold (all quantities are in dB): 

Moving from the model as described by equation (41). the hereafter derive and 
characterize statistically the observations that can be derived from level 
measurements in a mobile radio network and used for MS location purposes as used 
in embodiments of the present Invention. 

One measurement, or observation, as used In the embodiments of the invention for 
location purposes In algorithms A, B and C, is the difference between the average 
power received by the mobile, P'r, and the maximum power radiated by the i-th BTS, 

P t,max' 

= -^r - Pi,max (42) 
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This definition is justified by the fact that the MS location-dependent information Is 
embedded In the difference P'rP'i.inax and that the maximum radiated power P't,max is 
a parameter reasonably easy to calculate, given the network configuration. 

5 According to the propagation model (41), U(d'(x,y),v|/'(x,y)) is a random variable, due 
to the stochastic nature of u'(x,y). Since u'(x,y) is a Gaussian random variable, 
L'(d'(x,y),v|;*(x,y)) is a Gaussian random variable as well: 

L'(x,y)^J^(fii{x,y),ai(x,y)) ^43^ 

10 The mean value jJ.'L(x,y) and standard deviation a Ux.y) can be derived as follows: 
IVIean Value of L'(x,y) 

tJ'%(x,y) = ^[L'(x,y)] 

= -PV{d*ix,y))-APli^/;^(x,y)) 



Variance of L'(x,y) 

(4(a^,y))' = E [(L»(a;,y))'] - {A{x,y)y = E [{u'ix^y))'] = (al)^ 



15 



(45) 



The mean value and standard deviations of the i-th observation depends on the MS 
coordinates (x,y). The probability density function (pdf) of L' conditioned by x and y 
can be then expressed as 



20 



JL^x,y \x,y\ — y rexp < — -to— > 



(46) 



In embodiments of the present invention for MS location estimation purposes, level 
observations from different BTS's are used. For this reason, the covariance of two 
level observations is of interest. The covariance of two level observation in tum 
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10 



depends on the cross-correlation of the slow fading processes affecting the 
propagation of a signal from different BTS's. Often the slow fading processes are 
considered unconrelated. However, as Is known a certain conrelation exists between 
signals sent from different BTS's. Therefore in embodiments of the present Invention 
use a model for the cross-con-elation of the slow fading, defined as follows: 

(47) 

The Covariance of L'(x,y) and U(x,y) can therefore be shown as the following (x and y 
are neglected for notational convenience). 



Level measurements from multiple BTS's, BTS^ BTS"^ can be collected and used 

in further embodiments of the invention to estimate the MS location. Single level 
observations as shown In equation (42) can be stacked in a Nx1 vector of 
15 observations, L: 



L=[L\...,L^f 



(49) 



20 



which has a multivariate Gaussian distribution: 

The mean value and covariance matrix of L can in embodiments of the present 
Invention be readily calculated by the use of the results from the previous equations: 
The mean value of L(x,y) can be written as: 

mi,(x,y) = [fjiix,y), . . . {x,y)f ^ 



25 



wo 2004/023155 



PCT/IB2002/003644 



25 



10 



15 



The covariance matrix of L(x,y) can be written as: 

Ri,(x,y) = E {lL^} - mirni.^ ^g^) 

The generic element of RL(x,y) being (see equation (48)) 

^ { 4^W(a:,y) Mi (53) 

The probability density function of L conditioned by x and y can therefore be 
written as equation (54) (where |RL(x,y)| indicates the determinant of the 
covariance matrix RL(x,y)) 

fn^,y (L|a:,y) = '^~ym'f^^^^~^ {~\ ^ ~ ^T^i^>y)f'^^~\^>y) [L " mL(a;,y)]J 

(54) 

In further embodiments of the invention, one simplifying assumption can be 
made by considering the slow fading from different BTS's as having the same 
variance: 

i j A 



This assumption ' is accurate where the propagation takes place in a 
homogeneous communications environment. In other words where the 
communications environment is consistent and similar Thus the slow fading 
affecting different BTS's have statistically the same properties. The assumption 
20 of the same covariance for all slow fading links results in a modified structure 
of the covariance matrix RL(x,y). The structure of the covariance matrix 
becomes the product between the slow fading variance (which is common to all 
BTS's), au^, and the matrix, rL(x,y), which depends only on the location- 
dependent cross correlations, pu^'^(x,y): 

25 (56) 
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Using the assumption defined above, tfie probability density function of L 
conditioned by x and y in equation (54) can be written as: 

^ fM..y (L|a:.y) = (2^)^/2^ j^^^^^^), 1/2 {-^ " mx.(a:,y)]^rx.-(a:,y) [L - mi.(a;,y)]} 

(58) 

In furtlier embodiments, the slow fading processes are assumed not only to 
have equal variances but to also be uncorrelated. In such an embodiment 
10 rL(x,y) becomes an NxN Identity matrix, Independent of the MS coordinates: 

rL(a:,y) = Tl = I (59) 

and the probability density function (58) changes In 

15 

where (see equations (51 ) and (44)) 

1|L - mL(a;,y)|p = [L - mi.(rc,y)]^ [L - iai,{x,y)] = E + + APl{x,y)f 

(61) 

20 A l\/lodel for the Level Difference Observation 

The l-th level observation U defined in equation (42) Is obtained from the level 
measured by the MS on the signal transmitted by the l-th BTS. In a further 
embodiment of the present Invention observations can also be derived from the 
difference between the level measured by the mobile device from one BTS talcen 

25 as reference and the level measured by the mobile device from another BTS. 
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Using difference in levels help in eliminating unknown common biases in the 
absolute level measurements. Following an approach similar to those discussed 
above, the following shows the statistical characterization of the level difference 
observation. 

5 

The difference in level observation is defined in a further embodiment of the 
present invention as the difference between the level observation from the i-th 
BTS and the level observation from a reference BTS (in the following identified by 
the index i = 1) as defined in (42): 
10 D^ = L^-P ; i = 2,...,N (62) 

D' is a Gaussian variable, due to the stochastic nature of the log-normal fading: 

D' (x, y) ~ N(ni, (x. yl aj, (x, y)) (63) 

15 The mean value and covariance of D'(x,y) can be derived as follows: 
The mean value of D^(x,y) 

fJ-hM = E[D*{x,y)] = E[L\x,y)-P{x,y)] 

= -[PL'{d'{x,y))-PUid'{x,y))] 

- [API mx,y)) - API (64) 

20 The covariance of D^(x,y) and D^(x,y) 

The covariance between a pair of level difference observations is needed In 
embodiments of the present invention because the location algorithms use joint 
multiple level difference observations. The general definition of covariance 
25 between D'(x,y) and D^(x,y) is (x and y are neglected to simplify the notation): 
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The term resulting when 1 = j is the variance of the l-th difference observation 
(o'd)': 



= E 



{L'f] + E [(LO'] - 2E [L^L^] - 
5 = K-O +2criai(l-pM) ^^^^ 

Where p^'u is the cross-correlation between the slow fading affecting the 
propagation from BTS* and BTS' defined In equation (47). 
The term for i^t j can be calculated as follows: 

nD'Di]-t,^^,J^ = E[(L^-L0(L^-Li)]-(;*i-^i)(^i-4) 
= E [(L^)"] - E [L^I}] - E [L^U] + E [L'U] + 

+ {<^ipi^ + u^li^l) - {lA - /^D {A - A) 

(67) 

Thus in embodiments of the present invention the covariance between two level 
difference observations is summarized as: 



(68) 



The probability density function of d' conditioned by x and y can be expressed as 
follows: 



^ ^ V^c7l,(a:,y) ^| 2((7i,(x,y))" / 



(69) 
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Measurements from multiple BTS's, BTS^..•,BTS^ can be collected and used to 
estimate the MS location. In further embodiments of the present Invention Level 
difference observations can be stacked in a (N-1)x1 vector of observations, D: 

5 ^ ' ' J (70) 

which has a multivariate Gaussian distribution: 

D ^ ^(mi>(x,y),RD(x,y)) (7-j ) 

10 Mean value and covariance matrix of D can be readily calculated from the results 
defined above; 
The mean value of D(x,y) 



15 The covariance matrix of D(x,y) 

(73) 



T^i>{x,y) = E {DD^I - mDiHD^ 



The generic element of RD(x,y) being (see equation (68)): 

20 

The following probability density function of D conditioned by x and y results 

Mx,y (D|a:,y) = ^^^^^ ^^|i/2 ^ {-\ " ^t>Mf 'R'D~^{x,y) [D - mi>(a:,y)] J 

(75) 
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In further embodiments of the present invention the slow fading from different BTS's 
are assumed to have the same variance: 

K = cri = <Tu ; iJ = 2,...,iV (76) 



The elements of the covariance matrix RD(x,y) become: 

[R^(xv)] -/ 2a„2(i-pM) ; i 



= J 
¥^3 



(77) 



10 



Moreover, in further embodiments of the invention if the slow fading processes are 
assumed uncon-eiated. In equation (77) pv,^-' = po^J because /J^tl and po'"' = 0 when 
i^j; thus the generic element of Rd becomes: 



(78) 



The covariance matrix can be written as 

2 1 ••• 1 



1 2 
1 ••• 



••• 1 

1 2 



= a«2{l + i} = a„2ri5 



15 



where I and 1 are two (N-1)x(N-1) matrices 

• 0 



1 = 



1 0 
0 1 



■•• 0 

0 1 



; 1 = 



1 • 

• • 

1 • 



(79) 



(80) 



Using the assumption of slow fading uncorrelated and with equal variance, the 
20 probability density function (pdf) (75) can be written as follows 
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/d|x.v (P\x,y) = (3^)iv/2^^)^^ji/2 exp {-— ^ [D - xn^ix,y)f r^-' [D - m^ix,y)]] 



(81) 



The determinant of matrix m and its inverse can be calculated explicitly: 

(N-l) -1 ... -1 
-1 (N-l) ••. : 



-1 

-1 



= 1-^1 



(82) 



Substituting definitions (70) and (72) and the expression of ro'^ just obtained in the 
argument of the exponent in (81) it results: 



[D-mD(a:,y)f rD-MD-mD(x,i,)] = [D" - /Ji%, . . . - fi^^ [l - ^] 



2 



E - Mj>(x,y))' - {^D^ - A4,(x.y)^ 



10 



(83) 



The location algorithms used in embodiments of the invention are applications 
of the Maximum Likelihood (ML) principle to the level and level difference 
observations defined above. The Maximum Likelihood principle as discussed 
15 previously is a method widely used in the estimation theory. A brief review of 
the maximum likelihood principle follows. 

If X sD is an unknown random parameter defined in a certain domain D and y is 
an observed random parameter. The principle of the Maximum Likelihood (ML) 



wo 2004/023155 PCT/IB2002/003644 



32 

provides an estimate of x by maximizing the joint probability density function 
(pdf) of X and y: 

xjjfi = argmax/(x,y) 

x€0 (84) 

The ML estimates can be also calculated by maximizing the natural logarithm 
of f (x.y), A(x)= In f (x.y); usually referred to as the log-likelihood function: 

xml = argmaxA(x) 

*xei> (85) 

By writing f(x,y) as the product of the a posteriori probability density function of 
the observation y given the unl<nown x, f (ylx), and the prior probability density 
function for the unknown x, f(x). In other words the probability density function of a 
first unknown x, f(x), multiplied by the probability density function of the observation y 
conditional to the first unknown x, f(y|x). A similar expansion of the log-likelihood 
function becomes A(x) = In f(ylx)+ln f(x). In both cases the ML estimates can be 
calculated in embodiments of the present invention by solving the following 
problem: 



XML = argrc^ {In / (y I x) -I- In / (x)} 



(86) 



If no prior probability density function for the unknown x is available, f(x) can be 
neglected (or equivalently, x can be assumed to be uniformly distributed over 
the domain D), resulting in the simplified form of equation (86) shown in 
equation (87). 

XML = arg max {In / (y | x)} 

""^^ (87) 



In the context of MS location with signal level measurements, the observation y 
can be either the vector of level observations L defined in (49) 
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or the vector of level difference observations D defined in (70): 



(89) 



The parameter to be estimated x comprises in embodiments of the present 
invention the estimated MS coordinates (x,y) along with the estimated variance 
of the slow fading cu^: 



(90) 



It is shown above that both level observations L and level difference 
observations D can be modelled as multivariate Gaussian random variables. In 
such circumstances the ML criterion (84) can be further simplified. If the 
observation y eSl**" conditioned by the unknown parameter x is a M-variate 
Gaussian random variable with mean value my eSR*^'" and correlation matrix Ry 
gg^MxM_ probability density function of observation y conditioned by x is 

^^'""^ ^ (27r)^/^|R^(x)|^/^ ^""P {~l ~ ^y(^)f^-'i^) [y - my(x)]} 

(91) 

where the dependence of Ry and my from x is explicitly Indicated. The natural 
logarithm of fyix(ylx) is 

In/ (y I x) = -:^lu27r - |ln|R,(x)| - 1 [y _ my(x)f R^-^x) [y - my(x)] 

(92) 

and the ML estimate of x can be calculated according to the ML criterion (89) as 
follows: 
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XML = argmm {in |Ry(x)| + [y- my(x)f Ry-i(x) [y - my(x)]} 

(93) 

Equation (93) provides thus ttie ML criterion for MS location estimation when 
level observations and level difference observations comply with the Gaussian 
statistical models derived above. For simplicity, the following example Is 
restricted to the cases In which the variance of slow fading from different BTS's 
is assumed equal (the statistical models for the observations are described 
above). 

When N level observations are considered and the variance of the slow fading is 
assumed equal for all observations, the ML criterion for a Gaussian unknown 
(93) applies with the following definitions for y. my(x), Ry(x) and M in 
embodiments of the invention: 



where rL(x,y)=— Ru(x,y) and RL(x,y) is the covariance matrix of the level 
observations L. The ML criterion to estimate x, y and is thus: 



my(x) i- m.i,(cc,y) 
Ry(x) i- (Xu^rL{x,y) 
M = N 



y ^ L=[l\...,L^Y' 



(94) 




X 



(95) 



When N level observations are considered and the slow fading is assumed 
uncorrelated with equal variance for all observations, the ML criterion for a 
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Gaussian unknown (93) applies with tiie following definitions for y, my(x), Ry(x) 
and M: in embodiments of the present invention. 

y ^ L=[L\...,L^f 

my(x) -f- ntLix,y) 

Ry(x) (Ttt^I 



5 The IVIL criterion to estimate x, y and becomes: 



X 

y 



= arg 



min |lncr«2H-J:^||L-mi.(a:,y)f} 



X 

y 



(97) 



10 



15 



or analogously by Inserting the values of ||L-mL|| given in equation 61. 



X 

L y 



= ajg min 

X 



iin LaJ + J^^ {U + -pui^^y) + APi{x,y)f\ 



(98) 



In embodiments of the present invention <s^^ can be estimated separately. For 
fixed X and y this estimation of au^ is equivalent to finding the (strictly positive) 

minimum of the function f(s)=lns+— , which is Smin=K. The estimated value of 



au^ results in 



= E + PL'(^>J>) + AFl{x,y)f 

«=1 



(99) 



The ML estimation of x and y can be found in this embodiment of the invention 
by solving the following minimization problem: 
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= arg man 
X 1 




(100) 



where Dxy is the domain of existence of x and y. 

The domain Dxy can be determined In further embodiments of the present 
5 invention, for example, by using any of the additional location information 
available for implementation of Network Based Software Solutions (NBSS's) For 
instance, if the Cell Identity of the serving BTS is known, Dxy can be defined in an 
embodiment of the Invention as the geographical region served by such BTS. If In 
addition also the Timing Advance is available, Dxy can be defined in a further 

10 embodiment as the geographical region determined by the intersection of the serving 
area of the BTS identified by the given CI. In a further embodiment Dxy can be defined 
as a circular crown or ring with origin at the serving BTS coordinates and Inner/outer 
radii detemnined according to the TA value, for example by using the techniques 
known in the art. Restricting the domain of x and y in embodiments of the Invention, for 

15 example according to CI and TA, has two advantages: the first one Is that the TA 
infonnation is implicitly taken Into account; the second one is that the convergence of 
the minimization algorithm is made faster. 

When N-1 level difference observations are considered and the slow fading is 
assumed uncorrelated with equal variance for all observations, the ML criterion for 
20 Gaussian unknown (93) applies with the following definitions for y, my (x), Ry (x) and M: 



my(x) 

Ry(x) 

M = N-1 



■f- iaD{x,y) 



(101) 



where ro is the following (N-1)x(N-1) matrix (independent on the MS coordinates): 
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ro = 



2 1 
1 2 



■•. 1 
1 2 



(jv-i)x(^r-i) 



(102) 



The ML criterion to estimate x, y and ctu^ can be derived by using the result (83): 



X 



= arg mm < 



i=2 \j^2 J 



N 



i=2 



(103) 



where 



and 



ly^L"-- U (j = 2, . . . ,iV) 



(104) 



10 



15 



Mi>(^,y) = - [PL^ (d\x,y)) - VD (d^(a:,y))] - {API i^\x,y)) - API {^ix,y))] 

■T 

(105) 

Dxy is the domain of existence of x and y. Several possible definitions for Dxy 
are given later. 

Restricting the domain of x and y, for example according to the values CI and 
TA, has tv^o advantages: the first one Is that the TA information is implicitly 
taken into account in the estimation; the second one Is that the convergence of 
the minimization algorithm is made faster. 



20 
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In further embodiments of the present invention alternative interpretations of the 
IVIL criterion exist. f(x,z) the joint probability density function of 
5 unl<nown x and observation z can be written as product of the conditional 
probability density function f(xlz) and the marginal probability density function of 
the observation f(z). The log-likelihood function can thus be written as A(x) = In 
f{xlz) + In f(z) and the ML estimates is calculated by solving the following 
problem: 

^ML = argm^{ln/ (x | z) + In/ (z)} = argmaK{ln/ (x | z)} 
10 (106) 

The second equality holds because f (z) does not depend on x. If the Maximum 
Likelihood estimate of x is Inside the domain D, then it can be calculated as being 
the root of the following equation: 

^ln/(x|z) = 0 
15 Sx •' ^ ' ^ (107) 

The location method used in further embodiments of the present invention are 
based on maximum likelihood criterion as defined in (106) where the observation 
z comprises CI and TA information from the serving BTS and Cl's and RXLEV's 
20 (received level) values from all BTS's involved, and x comprises the unknown 
coordinates of the MS: 

x=[x,yf ^^08) 

In location serviceapplications, the knowledge of the CI of a certain cell implies that the 
25 geographical coordinates of the BTS antenna, as well as other parameters such as 
antenna orientation, cell width, transmitted power, etc, are known. 
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The ML estimate of the MS coordinates is determined in further embodiments of 
the present invention by solving the following minimization problem: 

{x,y) = axg max In / (a;,y | z) 

5 or, alternatively, by calculating the roots inside the domain D of the following 
simultaneous equations: 



To apply the ML principle, in the form expressed above, the domain of the solution 
10 D and the conditional probability density function f(x.ylz) need to be determined. 
These are detailed below. 

D is the domain where the location methods described in this document look for 
the solution x = [x,y]'^. D in (109) can be defined by using some a priori 
15 information on the region where the handset is possibly located- Several 
possibilities exist; four such methods used in embodiments of the invention are 
described below. 

1. D determined from Cell Identity (CI) of serving cell, 

20 In this case, D represents the geographical region where handsets connected to 
the serving cell are most likely located. D can be thus defined as the confidence 
region associated to a location estimate based on the serving cell CI information. 
A method to determine such a confidence region is detailed later. Figure 8 shows 
such a confidence region. The confidence region comprises the difference in area 

25 between two, unequal radii, circle segments 801, 803. The two segments 801, 803 




(110) 
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with a common origin 805 and common arc angles 807 are defined by tlie 
following set of parameters : an origin 805 located at point with coordinates 
(Xo.Yo), an inner radius Ri, an uncertainty radius R2, a orientation angle 809 a 
and Inclusion angle 807 p . The orientation angle 809 defines the angle from the x 
5 axis to the start of the arc. The inclusion angle 807 defines the angle from the start of 
the arc to the end of the arc. 

2. D determined from Cell Identity (CI) and Timing Advance (TA) of serving cell. 
A method to determine the parameters of the confidence region as shown in figure 
10 8 when, both the CI, and the TA information from the serving BTS is available (the 
TA, in particular, affects the radii Ri and R2). When the TA from serving cell is 
available, D can be thus determined with the method provided later eventually 
neglecting the cell sectorization (i.e., assuming a = 0 and p = 2n). 

15 3. D determined from Cell Identity (CI) of all cells involved in location calculation. 
The coordinates of the BTS's involved in the location estimation provide 
themselves an indication of the geographical region where the MS is located. D 
can be thus defined, for instance, by the convex polygon having vertices at the 
coordinates of the outermost BTSs involved In the location calculation. The 

20 concept is shown in figure 9. Figure 9 shows six base transceiver stations (BTS) 
with the boundary of D defined by four BTS's 901, 903, 905, 907. and therefore 
defining a quadrilateral 909 and two BTS's 911, 913. located within this D region 
boundary 909. 

25 4. D determined from coverage prediction maps. 

Following the same criteria proposed in application PCT/EP0 1/01 147 to determine 
the confidence region of a location estimate, coverage prediction maps for serving 
and/or neighbouring cells can be used to determine D. If the TA information from 
the serving BTS is available, the circular crown centred at the serving BTS 
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coordinates with inner radius Ri and outer radius Ri + R2 (as shown in figure 8) 
can be used in addition to coverage maps in determining D. Figure 8 shows this 
area as the area difference between 

5 Expressions for the probability density function f(x,ylz) 

Several expressions for the probability density function f(x,ylz) are detailed 
below. In embodiments of the present invention the observation z comprises 
information derived from received level observations (RXLEV). 

10 Given the following definitions: 

• N is the number of received level observations (RXLEV's) used to estimate 
the MS coordinates; 

* P\ P'^r are the received level observations (RXLEV's) measured by the MS 

from the N BTS's involved, expressed in decibel- 

15 P\.max .P^.max are the maximum radiated power values from the N BTS's 

measured in decibels. The i-th maximum radiated power, P*t,max. represents the 
maximum power at the output of the i-th BTS antenna in the direction of maximum 
gain. The value P^max as defined in equation (4) comprises transmitted power, p't, 
maximum gain of the BTS transmit antenna, G^max, antenna losses, cable losses, 

20 etc.; all measured in dB.The total (positive) attenuation experienced by the signals 
transmitted by each of the N BTS's involved while propagating toward the MS can 
be expressed in decibels as 

z'^Plma.-P'r (i = 1, • • • ^^^^^ 

25 The vector z = [z\...,z^]^ represents the observation based on which the MS 
location is estimated according to criteria (109) or (110). The probability density 
function of interest here is thus 
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10 



By using Bayes Theorem. f(x,ylz) can be expressed as follows: 

/(.,y|z) = / (z|x,,) = ,.00 

L L f(-MfMdxdy ^^^^^ 

f(zlx,y) needed in (113) can be evaluated by using the results from the above 
section, where the so-called multiple level observation L = [L\...,L'^]^ 
corresponds to -z defined In the present document. The most general expression 
forf(zlx,y) = fL|x.y(-Llx,y), is 

^ ^''^'"'^^ ^ (27r)^/^|R,(a;,y)|V^ ^{"^ " ^^Mf^-^x.y) [z - m,(a:,y)]} 

(1 14) 



where R2(x,y)=E{2z"^} - mzmz^ = RL(x,y) is the covarlance matrix of z. mz(x.y) = 
E{z} =[PL\d^x.y)) + APV(M/kx,y))....,PLN(dN(x.y)) + AP\(v|/Vy))f = -mUx.y) is 
the mean value of z(x.y). PL'(d'(x,y)) and AP'tr(V(x,y)) are the combined Transmit- 
Receive antenna pattern and path-loss associated to the signal transmitted by the 

y'-y 

15 i-th BTS, respectively, \|/'(x,y) = tan"^ x*-x- 's the angle of arrival (AOA) of the 

same signal, d'(x,y)=-y/(x' -xf +{y -yf is the distance between the MS and the 1- 
th BTS and (x\y^) (x'^.y'^) are the x,y coordinates of the N BTS's involved. 

In absence of any other a priori information, the MS coordinates can be assumed 
20 uniformly distributed in D. This assumption leads to a joint probability density 
function of x and y defined as follows: 

f(x,y) = I ^ (^'2/) ^ ^ 
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where M(D) is the size of D, and D can be determined using one of the methods 
described previously. 

Substitution of (115) in (113) leads to the following expression for f(x,yl z): 



which can be evaluated by using (112). 

In a practical implementation it is very difficult to use (116) in the minimization 
problem (109). For this reason, alternative approximate definitions for f(x,ylz) 
need to be determined. f(x,ylz) can be written as product of the probability density 
functions conditioned by each single measured attenuation .f (x,ylz^ ),..., f('x,yiz'^): 



The i-th probability density function in (117) represents the lil<ellhood of (x,y) given 
the attenuation measured from the i-th cell, z\ Physically f(x,ylz') represents the 
spatial distribution ^of the IVIS when the signal received by the MS from the i-th 
BTS experiences an attenuation z'. One embodiment of the present invention 
described below, defines f(x,ylz') in such a way that the x,y coordinates are 
uniformly distributed throughout a region, and defined by the observation z^ and/or 
the radio coverage properties of the i-th cell. A further embodiment of the present 
invention also described later, does not assume the handsets are uniformly 
distributed but to instead assumes that the i-th cell is omni-directional, so that the 
handsets are assumed to be uniformly distributed angularly from the BTS site but 
not radially. 




ix,y) e V 



(116) 



N 



(117) 



f(x,y z'): uniform distribution over a region Di 
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When the MS is assumed to be uniformly distributed in a certain geographical 
region Di associated to the l-th cell involved, f(x,ylz') has the following expression 



where M(Di) is the size of the region Dj. Several possibilities exist to determine Di; 
1. Dj from coverage maps 

If coverage maps generated by coverage prediction tools are available, Di in 
(118) can be defined as the Hearability area of the i-th cell. Hi. The hearability 
area identifies the geographical region where the signals radiated by the i-th 
BTS reach the handsets with a signal strength that is above the MS sensitivity 
level. 

The definition (118) takes into account only the identity of the I-th BTS but does 
not use the actual observed attenuation. One further embodiment of the 
invention refines the definition of f(x,y|z'), by Including the measured attenuation, 
is to substitute the hearability area with the coverage area of the i-th cell, Ni, 
which identifies the geographical region where the signals radiated by the i-th 
BTS reach the MS with the attenuation observed from the i-th BTS. From a 
practical point of view, it could be beneficial to consider a range of values for 
attenuation Instead of a single value; for example z'+Az', to take into account 
with Az' of the attenuation's random fluctuations. In this case f (x,ylz^) has the 
same definition as in (118), with Hi substituted by Ni. 




(118) 



2. Analytical expressions for Di 

In absence of coverage maps, Di in (118) may be expressed with many 
analytical functions. One possibility is to define Di as a generalized sector as 
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Shown In figure 10. Figure 10 shows the area defined by Di as the area defined 
by the area defined by a segment of a circle 1001 with an origin given at 
coordinates of the i-th BTS (x'.y') 1003, a sector cell orientation 1005, ^\ a 
sector width of twice the angular width 1007 A<|)'. and a (front) radius 1011 R'f. 
The area defined is also that Included within the area of a smaller circle 1009 
with the same origin 1003 and a (back) radius 1013 r'b. The area Di can 
therefore be defined as follows: 

7>. . / - ^0'' + (y - y'? <R^ ; o < \i^'ix,y) - <l>'\ < A4>' 

where v|/'(x,y) = tan~' ^i~^ J's the angle of arrival of the signal transmitted by the 
I-th BTS. 



With the definition (119) for Di, M(Di) in (118) is equal to (R'F)^A<|»'-Kir-A(|»')(R'B)^. 
In its most general definition, Di represents a sector cell, but can also represent 
an omni-directional cell, by setting A<j)' = n and R'b = 0. 

In an alternate embodiment of the present invention, D| can be defined as an 
ellipse with centre at a reference point of coordinates (x'R,y'R) and with semi- 
axes ax,i and cry.i: 

Vi(x,y) : + Z^-' < 1 

'^x.i ^y,i (^20) 

With the definition (120) for Dj, f(x,ylz') In (118) Is a three-dimensional cylinder 
with constant height }^ = and elliptical base. The general elliptical 

M(Di) TKT^jCTy^ 

definition allows an approximation to the coverage of sector cells to be made, 
but can also be applied to omni-directional cells (by setting the origin (x'R,y'R) at 
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the l-th BTS coordinates and ctx.i = ay.i = r'f in (120), Di is Identical to tine one 
obtained with (119) for A<|>' = n and R'b = 0). 

f(x,y|z'): non uniform distribution for omni-directional cells 

5 

In further embodiments of the present invention the probability density function 
f(x,yl2') can be calculated from the joint probability density function of the distance 
between the MS and the i-th BTS, d', and the angle of arrival from the same BTS 
V(x,y), f(d',Vl2'). as follows: 

0 \ X -x\ J ^^21) 

Where 



15 



20 



25 



\JM\ = 









dx 


dy. 










dx 


dy 





^(x - a;*)2 + (y - yi)2 

(122) 



In the following example the case of the omnl-directional cell is considered. In this 
example the MS angular coordinate \|/' can be assumed to be Independent from 
the radial coordinate d' and uniformly distributed about the region [-n,%]. Using 
these assumptions the following simplified expression for f(d',\|/'|z') can be defined: 
/ (cf ,^'1^0 = ^/ ; m<n,d^>0 ^^23) 

Strictly speaking, the l-th observed attenuation z' depends in general on the 
distance between MS and the i-th BTS, but also on the gain of the transmit antenna 
in the direction of the MS. If the cells are omni-directional, as assumed here, the 
antenna gain contribution can be neglected, since BTS antenna radiate uniformly in 
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all directions and the attenuation can be considered as the function of the MS-BTS 
distance only. In this sense, the i-th attenuation is equal to the i-th so-called path- 
loss PL*(d'), dependent only on the distance between MS and the i-th BTS: 



Thus in summary, the probability density function f{x,ylz'), with omni-directional cells 
can be obtained by inserting (122), (123), and (124) in (121), as follows: 



^ ^ 27r \ V 1 y ^(a; _ a;.)2 + (y - y* 



(125) 



10 



15 



20 



It has been shown in experimental and theoretical studies that the average 
received signal power decreases logarithmically with the distance between 
transmitter and receiver, both in indoor and In outdoor environments; thus the 
path-loss at a distance d ^do in decibels can be expressed as 



PL(d) = PL(do) + 10nlogio(^)-l-u ; d>do 



(126) 



where n is the environment-dependent propagation exponent, do is the "close-in 
reference distance" and PL(do) is the average path-loss experienced at a distance 
do from the transmitter. In free space the value of n is 2 but the value n grows when 
the density of obstructions increases. Table 1 lists typical path-loss exponents in 
different environments. 



Environment 


Propagation exponent, n 


Free space 


2 


Urban area cellular radio 


2.7-^3.5 


Shadowed urban area cellular radio 


3-^-5 


In-building line of sight 


1,6-M,8 


Obstructed in-building 


4-^6 


Obstructed in-factories 


2-^3 



Table 1 
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do must also be selected on the basis of the environment. When cells are large, do 
is usually set to 1km; in case of micro-cells, the "close In reference distance" Is 
usually smaller (It can vary from 1m to 100m). PL(do) is therefore calculated In 
embodiments of the Invention using experimental data. When this Is not possible, 
5 PL(do) can be estimated in further embodiments by using the free-space path 
loss law; If do is close enough to the transmitter the idealized condition of 
propagation in free-space can be assumed {X = c/f is the signal wavelength, c is 

the speed of light and f the frequency): PLfree space(do) =101og 

10 In the model (126) u represents the shadow fading affecting the signal transmitted 
by the i-th cell. It is generally modelled as a random variable with log-normal 
distribution with standard deviation Ou (i.e., u measured in dB is a Gaussian 
random variable: u ^ N(0, au)). Typical values for Ou range from 5 to 10 decibel. 
By defining 

15 A = PL(do)-10nlog,o do .^27) 

B = 10n 

the model (126) can be re-written as. 

PL(d) = A + S logio d + u (1 28) 

20 (A well known model for the path-loss compliant with (128) is the Okumura-Hata 
model, where d Is the distance between MS and BTS measured in kilometers, A and 
B are functions of signal frequency, f, base station effective antenna height, he, 
mobile terminal antenna height, hm, and City Type (either "Large City" or 
"Small/Medium City"). The Okumura-Hata model is coherent with the formula (126) 

25 provided that do = 1km, A= PL(do), and B = lOn.) 




For a given path-loss PL, the distance of>do has the following expression 



wo 2004/023155 



PCT/IB2002/003644 



49 



d = 10 B 



(129) 



(130) 



For a given path-loss value and using tiie assumption of log-normal slow fading, 
tine probability density function of d in (129) can be calculated with known 
5 standard random variable transformation techniques: 

where C(do) is a normalization factor introduced so that •C/(^|P^)^/^= 1= 

C(rfo)=/ /(plPL)dp . 

J do (131) 

10 

Inserting (130) in (125). the probability density function of x and y for a given 
observed attenuation from the i-th omni-directional cell (see (124)) is as follows: 

/ I A _ B'/C'jd^) {-^ (g*logio - x'Y + (y - y^Y + Aff] 
V^'^'"' ) (27r)^/2 ai InlO {x - t^Y + (y - yi^ 

(132) 

j5 valid for d^{x,y) = ^/(o; - x^Y + (j/ - 2/^)2 > rfp. 

Figure 11 shows a series of plots of the probability density functions f(x,y|z') 
using equation (132) for different values of attenuation z' =PL'. Figure 11 shows 
that the probability density function has circular symmetry around the vertical axis 

20 going through the I-th BTS coordinates. It can be seen that, as the path-loss 
increases, the probability density function spreads.. At low values of attenuation 
the probability density function curves peak close to the l-th BTS and rapidly fall 
towards zero as the distance from the I-th BTS increases. At higher values of 
attenuation the peaks of the probability density function curves move further away 

25 from the I-th BTS coordinates. As the value of attenuation increases the 
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distribution becomes flatter - I.e. the peak probability density function value Is 
smaller but the rate of increase and decrease of the probability density function is 
lower. This spreading produces the suggestion that at higher levels of attenuation 
the probability of handsets Is further away from the BTS site grows. 



f (x.ylz'): empirical Gaussian distribution for omnl-directional cells 
In embodiments of the present Invention the probability density function f(x,y|z') 
(132) is approximated by a bl-variate Gaussian probability density function defined 
as follows: 



h = 2;^p^ exp |-i [x - mf R-^ [x - m]} , 



(133) 



where 



X = 



.1 

y\ 



m = E{x|:2'} = 







E{j/|z»} J = 





R = E{xx^} - mm"^ = 



(134) 



(135) 



(136) 



To derive the Gaussian probability density function in (133), m and R (more 
precisely, its determinant IRI and its Inverse R"^) must be detemilned. It can be 
noticed that the average values of x and y with probability density function f(x,ylz') 
defined In (132) are 



A*=r.i = E{a;|«*}=a;* ; /Uj,.^ = E{y|z»} = 



(137) 



thus 
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10 



m = 



(138) 



To obtain an expression for tiie correlation matrix R in (136), tlie expected values 
E{x2|z'}, E{y2|z'}, and E{xyi2'} need to be calculated. 

Tfie analytic expression of E{x^|2'} can be obtained as (the second equality in the 
following is obtained by solving the integral in polar coordinates: 

E{x^\z'} = f f x^f(x,y\Adxdy 

0 

B^/cUdo) r+'r o r 1 / . .x2l 



(27r)*'''°crilnlO J-iT ^do 



(27 



(2t 

= (a;')2:^+2i 



(139) 



where % and I, represent, respectively, the first and the third integrals above: 



(140) 
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The following definition for the complementary error function is used. 



(141) 



5 With a very similar derivation, it can be proved that E{y^|z*} has the following 
expression: 

Finally, the analytic expression of the correlation term E{xy|z'} can be obtained: 



10 



(a;* + />cos^)(y' + pan ^) 



' + pan 4>) ^ L 



(143) 



Using the results (139), (142), and (143) in the definition for R given in equation 
(136), the following R results 



R = 



2; + (x')^(2;-l) a:V(^-l) 1 



(144) 



15 



The detemiinant of R being |r| = (l.)^ +(x')^Ii(X -l)+(y')'li(^ -l) and its inverse being 

IRI L -^'y'i^ - 1) ^ + - 1) J ■ 



(145) 



With the above results it is possible to define the Gaussian probability density 
function as shown in equation (139) as follows: 
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10 



+ {Ax - x') - y\y - y"))^' ^^^^^ 

In further embodiments of the present invention the use of the path-loss model 
(128) is extended to distances below the close-in distance (i.e., for d<do) or, 
alternatively, the close-in distance is made tend to do^O the above functions 
remain well behaving and the probability density function tends to zero as d 
tends to zero. In the limit do^O, the term % tends to 1 (in fact, (0.5)erf{c(-oo)} = 
1). C'(do)-»1. and 



1 0 
0 1 



(rfo->0) 

(147) 



where 



^ ^ 1 /^InlO / i <7flnl0'\l 



(148) 



Inserting the above result for do^O into (133), the following approximated 
15 Gaussian probability density function results: 



(149) 



In figure 12, the behaviour of I^^ as a function of the path-loss/attenuation z' Is 
shown. Figure 12 shows a graph of ho of range 0 to 8 against path 
20 loss/attenuation z' of range 110 to 145 dB. The graph plot resembles an 
exponential type plot with a ho of slightly above 0 at z' =110 dB rising initially 
slowly but increasing rapidly as z' passes 140dB. 
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f(x,y|z'): non uniform distribution for omni-directlonal cells 



5 The location algorithms In further embodiments may be applied to non uniform 
probability density functions in omni-dlrectional ceils, By tai<ing logarithms of the 
terms from equation (117), the minimization problem can be rewritten as 

ix,y) = axg ma^^ g In / {x,y | z') 



10 or, alternatively, 



[||ln/(.,.|.0=O 



51) 



As described above, equation (138), defines the probability density function 
f(x,yiz') in the environment created by omnl-directional cells, based on the 
15 logarithmic path-loss model (128). By using this definition for the probability 
density function f(x,ylz'), the following expressions for the partial derivatives can 
be found 

^In/ {x,y I z") = F'{x,y) (x - x') 

i 

^lnf{x,y\z') = F\x,y){y-y') 
^ ^ ^ ^ (152) 



20 where 
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L 2<7<2inl0 J 

(153) 

The MS location can be thus estimated in further embodiments of the present 
invention by solving iterativeiy the following set of simultaneous nonlinear 
equations 



N 



=0 

1=1 ' 

(154) 



Equation 154 defines algorithm D which is a non linear algorithm. 

f(x,y|z'): empirical Gaussian distribution for omni-directional cells 
Earlier two Gaussian approximations to the probability density function ^(x.ylz') 
were defined. The first approximation provided was fG(x,ylz'), In equation (146). 
By using such expression, the following partial derivatives result: 

I ^ In / (:..y I = [-^(:, _ _ _ ,y(^ _ ^.-^ jj 

(155) 

Where the determinant of R is given by |R| = (li)' +(xO'li(^ -l)+(yO'l..6 -l). The 
IVIS location can be thus estimated in embodiments of the Invention by solving 
iterativeiy the following set of simultaneous nonlinear equations: 
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N 

E 

1=1 

N 

E 

i=l 



IRI' 



^{(a:0^x-a;V(y-y')}]=O 



(156) 



10 



15 



Equation 156 defines algoritfim E whicli Is a non linear, Iterative algorithm. The 
second approximation for the probability density function f(x,ylz') given in section 
is fcoCx.ylz') in equation (149). This probability density function is obtained as a 
limit of the first Gaussian approximation, fG(x,ylz'), in the limiting condition do-»0. 
By using such expression, the problem simplifies noticeably; In fact the following 
partial derivatives result: 



(157) 



where nx.i = x', Hyj = y" and Iio, defined in (148), depends on the i-th attenuation z'. 
With the result above, the MS location estimate can be calculated in embodiments 
of the invention in closed form as follows: 



E^ 

i=i Ao 



N yi 

N . 



(158) 



Equation 158 defines algorithm F which Is a linear algorithm with a closed form 
solution. 



Extensions of the closed-form method 
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57 

In the algorithm defined by equation (100) and used in embodiments of the 
present invention the estimated x (and respectively, y) MS coordinate is obtained 
as a weighted average of the x (and respectively, y) BTS coordinates, the signal 
of which is received by the MS. This algorithm can be extended by means of the 
following formulation: 

N N 



(159) 



where w\...,w'^ are suitable weights assigned to each one of the N BTS's 
involved. 



The weights in the algorithm (158) are calculated as the reciprocals of the terms 
Iio....,Ino. The weights being calculated in such a way that the weights decrease 
as the attenuation increases, as can be seen also from figure 11. It is possible to 
suggest that the signals transmitted by BTS's closer to the MS undergo lower 
15 attenuation than signals transmitted by BTS's located further away from the 
mobile device. This suggestion can be accounted for by imposing the rule that the 
weights in the generalized closed-form algorithm (159) are defined in such a way 
that, If the attenuation of the signal received from the i-th BTS is low, then w* Is 
high, and vice versa: 

20 z'/^ ^ ^^gQj 

Hereafter follows three empirical definitions for the weights w\...,w'^ according to 
criterion suggested above as used in further embodiments of the present 
Inventions. 



1. Definition of Weights 
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The first definition of ttie weigiits wliich reflects the above rules comprises 
defining the i-th weight as the reciprocal of the attenuation experienced by the 
signal transnnitted by the I-th BTS: 



This definition of weights is further enhanced in further embodiments of the 
present invention by the introduction of auxiliary variable parameters to be 
determined by using experimental measurements. 

2. Definition of Weights 

A further definition for the I-th weight used In embodiments of the present 
invention Is to use the inverse of the estimated distance, d* , between the MS 
and the i-th BTS, obtained from the i-th observed attenuation z': 



The level of the signal received by a MS (and the attenuation as well) does not 
only depend on the distance between MS and BTS, but also on the gain of the 
transmit antenna In the direction If the MS. However, if the N BTS's Involved in 
the location calculation are omni-directlonal, the contribution of the BTS 
antennas can be neglected and the attenuation can be considered to be a 
function of the MS-BTS distance only. For this approximation, the i-th attenuation 
z" is equal to the I-th path-loss PL'(d'), and is dependent only on the distance 
between MS and the i-th BTS d' =^[x-x'f +(y-yf : 




(162) 



z* = PL*(fr) 



(163) 
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Neglecting the slow-fading in model (128), the following expression for the 
weights as given in (162) produces the following: 

ti;* = — = 10 = 10 — — 5 

(164) 

5 A further alternative can be found when the weights defined in (164) are 
adjusted to fit experimental measurements. Two observations, described below, 
can be made in order to produce more accurate weights to be found: 

Since the MS coordinates are estimated using (159) with weights given by (162), 
10 the absolute values of the distances are not strictly needed. In the algorithm 



defined in equation (159) only the ratios 




contribute to the location 



estimate. This means that if the same relative error is made in estimating each of 
the distances used in the location calculation, the resulting location estimate is 
not affected. This automatic error cancellation factor means that going from 
15 outdoor to indoor or from open area in a city to narrow city canyon should not 
greatly affect the location accuracy. This error cancellation factor Improves the 
accuracy of the location estimate and minimizes the need to use a highly 
optimised and parameterised path-loss model In embodiments of the invention. 

20 For a location estimate calculation using the algorithm (159) and weights (164) 
only the functional behaviour of the distance as a function of the received signal 
strength is needed. In a simplified example, if the base station antenna height, 
central frequency, maximum radiated power are the same for each BTS used in 
the location calculation almost all the terms of the Okumura-Hata path loss 

25 model cancel out and the resulting weights to be used in the location calculation 
simplify so that the weights can be found by the equation (165): 
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3, Definition of Weights 

A further embodiment of the present invention determines experimentally a 
functional relation between the weights, w'. and attenuation, . It Is possible to 
collect a sufficient amount of received level (RXLEV) measurements, with 
corresponding GPS coordinates to be used to reference the exact MS locations, 
and from this Information solve the reverse location problem: If the MS 

coordinates are known In equation (159), the weights w^ w^ can be 

determined as a function of the attenuation by minimizing the error between 
estimated MS coordinates and exact MS coordinates. This approach fails where 
experimental measurements are not available. 

The location methods described above apply to a network made from omni- 
directional cells, where the BTS's transmit antennas radiate isotropically in all 
directions. In sector cells, the assumptions made in an all omni-directional cell 
environment hold with a reasonable level of approximation error only when the 
MS Is located in the main lobe region (MLR). In other words when the mobile 
device is in the region illuminated by the main beam of the BTS (directional) 
antenna. A sector cell equipped with a transmit antenna having Half Power Beam 
Width (HPBW) of 60° has a MLR which approximately extends angularly ±60 degrees 
about the direction towards which the BTS antenna is oriented. 

In further embodiments of the present invention, the algorithms defined in 
equations (154). (156), (158), and (159) are applied if the N BTS's involved In 
the location calculation, out of the total number of values received by the MS, 
are those that radiate the MS with the main beam of their transmit antennas. 
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Unfortunately it is not possible to determine whether the MS is in the MLR of a 
certain BTS from a single BTS received value because, in order to obtain this 
information, the location of the MS is needed. However, in embodiments of the 
present invention it is possible to determine with a low error rate if the MS is in 
5 the MLR of a certain sector cell when the MS measures and compares the 
received levels (RXLEVs) from different co-located sector cells. In this case, as 
the BTS antennas of co-located sector cells have different orientations, the MS 
is most likely located in the MLR of the cell the signal of which experiences the 
lowest attenuation among the co-located sector cells. 

10 

With the l^nowledge gathered by the method above the approximation found from 
the use of omni-directlonal cells is accurate (e.g., and embodiments of the present 
invention may apply the algorithms defined in equations (154), (156), (158), and 
(159)) if the N BTS's involved in the location calculation are selected out of the 
15 ones measured by the MS according to the following procedure: 

1. Select all omni-directional cells measured by the MS. 

2. Select all sector cells, where no co-located cells are measured by the MS at the 
same time. 

3. Select the sector cell with the lowest attenuation out of the co-located sector 
20 cells, where co-located cells are detected by the MS. 

If one or more of the discarded co-located cells produce a received level value 
attenuation that differs from the selected cell by an amount by a certain range 
(say 5 dB), select these co-located cells also and the average of the attenuation in 
25 embodiments of the present invention or in further embodiments of the present 
invention use the attenuation values of all of the selected co-located cells. 



CI and CI+TA based algorithms 
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These algorithms will be described In more detail hereinafter. These algorithms are 
linear and rely only on serving cell Information. They both deliver a location estimate 
and a confidence region for the location estimate. 



These algorithms are considered to act upon a mobile device being serviced by 
one generalised sector cell with a coverage or hearability area defined 
analytically or by practical measurements. The generalised cell is a simplification 
of the effect produced by the BTS antenna angular gain to both transmit and 
receive as shown by the practical antenna gain shown in figure 5. 

As explained previously and as shown in figure 10 the generalised cell is defined 
by a series of parameters comprising. (|> the sector orientation In degrees 
from the x-axis, the sector angular width. Rp the sector front radius, 
and Rb the sector bacl< radius. 

The parameters defining the generalised sector may be found from coverage maps 
or coverage prediction tools. The Input information required by the method 
comprises, the Cell Identity of the serving BTS (CI), the serving BTS coordinates; and 
the map of the serving cell. 

The steps performed to calculate the generalised cell parameters from a coverage 
map of the cell comprise: 

1. Selecting the coordinates of the serving BTS using the CI: Xs, ys. 

Given the CI of the serving cell, the x - y coordinates of the cell are identified the 

CI, (Xs,ys), and are retrieved from the database of BTS coordinates. 



2. Selecting the serving cell map: S 
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Analogously as In the previous step, the CI is used to select the serving map, S, 
corresponding to the current CI. S is typically determined by a network planning 
tool which uses a large set of information such as BTS configuration parameters, 
3-dlmensional terrain maps, tuned propagation models, etc. and, in order to 
5 determine the serving area, takes into account also the presence of other BTS's in 
the network. 

The serving map represents the geographical region where the ceil identified by 
CI is serving. In an arbitrary x-y Cartesian system, S can be obtained by dividing 
10 the region into elements of area {Ax)x(Ay) and representing each element by the 
coordinates of its centre: 

S'{xn,yn} I n = l,.,.,iV5 (i66) 

The coordinates (Xn,yn) represent the centre of the n-th pixel of area {Ax)x(Ay) 
15 where the cell is serving. 

An example of serving area is shown in figure 13. Figure 13 shows a base 
transceiver station (BTS) 1301, a main coverage area 1303 and minor coverage 
areas 1305. The base transceiver station 1301 lies within the main coverage area 
20 1303. Smaller coverage areas 1305 are positioned adjacent to but not touching 
the edges of the main coverage area 1303. All of the coverage areas are divided 
into a plurality of elements 1309 by a grid system 1307. The total of all of the 
coverage areas dependent on the BTS is known as the serving area for the BTS. 

25 3. Determination of the coordinates of the centre of mass of the serving map: 

Xmc»ymc 

The (x,y) coordinates of the serving cell's mass centre are formally defined as 
follows: 
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(167) 



where M(S) is the area of the serving map: 

M{S) = / dxdy 

JS 



(168) 
5 

By using the divided elements defined in equation (166) the coordinates of the 
mass centre can be calculated as follows: 



2 ^''s t\s 



10 4. It is much easier to calculate the desired estimates in a polar reference system 
(pi,0i) originated in the serving BTS of coordinates (Xs.Ys). In this polar coordinate 
system pi Is the distance of the point (xi.yi) from the BTS and Gi is the angle 
measured counter-clockwise from the x-axis : 



15 



(170) 



In the polar reference system the serving map is represented by a set of points: 

•Spolar ■• {Pnfin} \ n = 1, . . . ,Ns. (171) 



5. Determination of the main direction for the entire serving map: emc 
20 The main direction Gmc provides an indication of the serving sector's bearing, <|)s- 
This angle Is used as a reference direction for normalizing the angular coordinates 
of the cell's coverage area elements. Q^c can be approximated as the orientation 
of the cell's mass centre from the BTS coordinates. 
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In the polar reference system (p, 0) the mass centre of the serving area has 
coordinates (pmcQmc) such that: 

PMC = yjixs- XmcY + (ys - VmcY ; ^mc = tan"^ ^^^^j 

5 6. To further simplify the actual calculation of the estimates the new polar 
reference system is rotated in such way that zero and 2ti angles are as far as 
possible from the most important serving area. After this rotation all the needed 
estimates can be calculated using a known simple sorting algorithm. 

10 The rotation of the m-th pixel with angular coordinate Om, is defined as follows: 



^1 em-en 
\ em — eR 



X ^ J ^m-eR + 2n i£0m-eR<O 



where 

eMc > TT 



^ \ euo - TT if 



(174) 



15 The resulting S^'s are In the range 0 < 5m < 2%. The directions 5m's closest to the 
main direction Gmc are associated to 8^ = ti and the directions furthest away from 
9mc are associated to 6ni = 0 and 6m s I-k. 

Next step is to sort the points in order of increasing distance from the serving 
BTS: 

20 {A,<yi}£i = sort{pi}. ^^^^^ 

At this position the serving map is represented by an ordered set of points: 

«5polar : ; n = 1, . . . ,JV5. 

(176; 



25 7. Determination of the front radius, Rp: 
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where [•"] rounds • to the nearest integer towards plus infinity (using the (<nown 
ceiling operation where all numbers are rounded up towards the next positive 
integer) and ,y' determines the desired fraction of points in set Spoiarto have the 
5 distance from serving BTS less than Rf. A good value for the parameter y is 
between 0.95 and 0.98, in some embodiments of the invention. 

8. Determination of back radius (Rb), orientation as a function of distance from 
serving BTS (<j)s(d)), and angular width as a function of distance from serving 
10 BTS (A(j)s(d)). 

First it Is helpful to define the parameters of the circular crown depicted in figure 
14. Figure 14 shows a constant timing advance (TA) region defined by the 
difference in area between two concentric circles with a common origin 1401. 
15 The first circle 1403 having a radius Ri 1407 and the second circle 1405 having 
a radius R1+R2. 

The mathematical definition for the crown in the polar reference system Is: 



Where Rinf is inner radius and Rsup is outer radius. In the example shown in figure 14 
the first circle radius 1407 Ri = R|nf., and the second circle radius Ri+R2=Rsup. 

This crown is therefore represented by a set of points: 



C:{{p,9)eIR'':B^nf<P 




(178) 



20 



25 



xrown • 



{pn,Sn} ; n = 1, . . . ,Ncromx{Ni 



crown 



<Ns). 



(179) 



9. To determine the back radius (Rb) the following steps are followed: 
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(a) set i = 0 

(b) The number of points in the crown defined in equation (179) is caicuiated using 
Rinf = i-A and Rsup= (i+2) A, where A is the distance of digitized elements in x 
and y directions. 

5 (c) If the serving area is omni-directional at the distance Rsup, there are 
approximately Njdeai = 7r(R%up - R^inf)/A^ points In the crown. If the actual number 
of points in Scrown. Ncrown is less than y'Nideai, , the serving area is no longer 
considered to be omni-directional. The value of y' < 1.0 and is based on 
experience a good value being 0.75. If the serving area is not considered to be 
10 omni-directional the back radius Is chosen to be Rb = Rsup- 

(d) If the serving area can be considered omni-directional at distance Rsup. set 
1=1+1 and go to (b), otherwise the estimate of Back radius is ready. 



10. To determine the orientation, <|)s(d), the angles in the circular crown are sorted 
15 in order of either increasing or decreasing angle; 



From the definition of the crown Ri„f = d-2A and Rsup = d+2 A. The estimate for 
the (|)s(d) is therefore the median value of the angles and it is obtained as: 

20 Md) = Si\i=n^„^/2V ^^g^^ 



11. The Angular width as a function of distance from serving BTS, A^s{d) is 
calculated using the same sorted set of points In the circular crown as is used for 
(|>s(d). Angular width. A(|>s(d), Is calculated as follows: 
25 ^Md) ^ max{\\6i - MMMSn^^ - MM)- (182) 
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12. As the last step it is necessary to rotate <|)s(d) and A<|)s(d) bacl< to the original 
co-ordinate reference system. This is done by simply adding 0r to the obtained 
values. 

5 Where coverage maps are not available, it is possible to analytically estimate the 
cell front radius, Rp, and cell back radius, Rb. with the method presented below. 

As described earlier, such as by equation (33), the average received power, Pr, 
by a mobile device a distance, d, from the serving BTS can be expressed as: 
10 PR{d) = Pt + G- PL{d) (1 

The BTS transmits a l<nown transmitted power Pt and the gain of the antenna 
installed at the serving BTS in the direction of the MS is G. PL(d) is the path- 
loss affecting the signal power as it propagates from serving BTS to MS. The 
15 accepted model, such as those used in equations (126) to (129) above, assumes the 
path-loss increases logarithmically with the distance d : 

PL{d)=A + B\ogd (184) 

The sector front radius can be defined as the distance at which the average power 
20 received by a MS, Pr, is above a certain threshold defining the cell edges, P\. The 
effect of potential shadow fading is allowed for by including a shadow fade margin 
FMff (ail quantities are in logarithmic units) 

25 Under the assumption of log-normal slow fading, the fade margin can be defined 
as FMa=za, where a is the standard deviation of the slow fading and z is such 

that F(z) ="/-Q(^z; = l~^|'e-'^'''dx is the radius estimate's reliability. 



wo 2004/023155 



PCT/IB2002/003644 



69 



According to equations (183) and (185) tlie front cell radius In kilometres can be 
then expressed as 



Rf = 10 B ^^3QJ 

wliere Gm is the maximum antenna gain in dB. 

The sector back radius can be calculated from the front radius by including 
information on the radiation pattern of the antenna Installed at the serving BTS 
site. Figure 5 shows an annotated version of figure 4, a plot of antenna gain in 
decibels against angle measured In degrees. The antenna gain, symmetrical 
about the angle G=0, comprises a main lobe 501, and four side lobes 503, 
505,507, 509. The main lobe 501 has a maximum value of Gm at 9=0, whilst 
none of the side lobes 503, 505, 507, 509 have gains greater than pG^. p is the 
maximum Back-to-Front Ratio, that Is the maximum ratio between the antenna 
gain In the direction of maximum radiation, G^ and the gain of the antenna in the 
directions outside of the main lobe, it Is evident that the cell front radius, Rf, is 
related to the maximum gain in the main lobe directions, Gm. and the cell back 
radius, Rb, is related to the maximum gain in the directions outside of the main 
lobe pGm ■ 

By using the simple propagation model (183) and the following rough 
approximation of the antenna gain G(e) in dB (Gm = 10iog(gm)) 



G{e) - * 



m 



0<e< OadB 

^3dB < ^ < 27r - ^3dB 

27r — ^3dB <9 <2ir (187) 



it Is possible to write the following upper limit for Rb: 

Rb< P^Rf (188) 
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Typical valjjes used in GSiVIQOO are Pr = 50dBm, P"'r = -95dBm. For a BTS 
lieiglit of 30 meters, typical values for variables are A = 124.5 and B = 35.7. In 
the mobile radio environment, a is often assumed equal to 8 dB. For this 
example the BTS antenna has HPBW 65 degrees, Back-to-Front Ratio -18dB (p 
5 = 0.0158) and maximum gain 12 dB. With these values and for z = 0.675 (i.e., 
F(z) = 0.75 equivalent to 75% cell radius reliability) formulas (186) and (188) 
give Rf = 5.7425 kilometers and Rb = 90 meters. 

From the examples shown, it is evident that the representation of the sector as 
10 shown in figure 10 might approximate poorly the real cell coverage. In particular, 
A(|)s is usually too large to represent the sector width in those regions far away 
from the BTS coordinates. To improve this inaccuracy in estimation, it is possible 
to define the serving sector width as a function of the distance from the serving 
BTS, A<|>s(d). By taking this into account it is possible to neglect In the location 
15 algorithm those regions which clearly lie outside of the serving area. The same 
applies for the sector direction, ^s, which can be assumed to be a function of the 
distance from the BTS site, (|)s(d). 

Rf. Rb, and (t)s(d) (eventually dependent on the distance, d) create a 
20 analytical estimation of the simplified borders of the serving cell as follows: 



In embodiments of the present invention the above estimations of the 
generalised cell parameters are used to estimate the location of the mobile 
25 device (MS). The algorithm used in embodiments of the present invention is 
referred to as a Classic CI-TA Location Algorithm, as opposed to a Map-Aided 
CI-TA Location Algorithm. 



= Rf ; 0 < \ip{x,y) - 4>s\ < A4>s 
= Rb ; |^(a:,s/) - </>s\ > A(f>s 



(189) 
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As explained previously the input parameters to tlie algorithm used in 
embodiments of the present invention are CI, TA (if available) and confidence 
coefficient of the location estimate, ^. As also explained previously the output of 
the calculation is a location estimate and a confidence region (e.g., a 
5 geographical region within which the real MS location is situated within a degree 
of confidence The result of the location calculation differs if the cell is 
sectorized or omni-directional. 

Embodiments of the present invention described below are capable of providing 
10 these results when both CI and TA are available, further embodiments also 
described below provide results when only CI is available. 

If the serving cell is sectorized, the location estimate is calculated by combining 
the serving BTS coordinates, Xs and ys, with an estimate of the i\/IS-to-serving 
15 BTS distance, d, and an estimate of the angular coordinate of the MS from the 
serving BTS site, vj>: 



d and \j/ are the estimated distance and angle respectively as described below. 

20 

Given the TA information, the method described previously can be used to 
determine an estimate of the distance between MS and serving BTS and the 
radii of a circular crown (C) centred at the BTS coordinates where the mobile 
station can be located with a confidence y- 



The distance estimate, d, is calculated as the 50-th percentile (or median value) 
of the real distance, d, i.e. d is such that 




25 
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Pr(d < d) = 1/2 



(191) 



Other analogous definitions for d , such as the mean of d, are possible. Using the 
results described in European Patent Application number 102251 which is hereby 
5 incorporated by reference, the median estimated distance is: 



lb = 3.69ias Is the bit period and c = 3 x 1 0° m/s Is the speed of light. 

In absence of any other information, the angular coordinate of the MS can be 
estimated with the orientation of the sector, I.e. vpr can be set equal to <|)s. 
15 If the sector orientation is a function of the distance from the serving BTS, then \j/ is 
the sector orientation at a distance equal to the one estimated on the basis of the 
TA, I.e. xfr = <|)s (d). 

Further means for estimating the MS angular coordinate can be Incorporated into 
20 the algorithm found in embodiments of the present invention; for example, y can In 
further embodiments of the invention be determined by processing signal level 
measurements (RXLEVs) performed by the MS or, in the future, by using the angle- 
of-arrival information made available by (smart) antenna arrays installed at the 
serving BTS site. 



Associated to the estimated MS location is the confidence region R as shown In 
figure 8 and comprising the following parameters: origin located at point with 



rf = <*rA + TAc ; TAc= -Xi/2 



(192) 
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coordinates (xo.yo); inner radius. uncertainty radius, Rz; orientation angle 
measured counter-clocl<wise from x axis, a, and Inclusion angle defining the width 
of the sector, p. 

Analogously to the distance estimate, the confidence region can be determined by 
using the method described in European Patent Application number 102251 which 
is hereby incorporated by reference. 

The origin of the confidence region is at the serving BTS site: 



Given the measured TA and the confidence coefficient for the distance estimate y, 
the confidence interval of the distance estimate (d) can be determined from the 
statistical properties of the TA measurement provided by the l<nown map aided Cl- 
TA estimation technique. The confidence interval for the distance estimate is 
determined by n and rg defined as 



where X(i±yy2 Is the 100—^% percentile of the TA measurement error, n and rg 

are such that the true MS-serving BTS distance, d, falls within the confidence 
Interval [d-n, d+rg] with a probability ,y (as shown in figure 15): 




(194) 




(195) 



Pr(f? - Ti < d < rf H- Ts) = 7 



(196) 
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Figure 15 shows a base transceiver station (BTS) 1501, a location estimate d 
1507 which lies within the area defined by two arcs 1503, 1505. The larger arc 
has a radius dsup, the smaller arc has a radius djnf. 

In embodiments of the present invention the confidence region parameters (Ri, R2, 
a, and p) are calculated slightly differently depending if the serving cell back radius 
(Rb) is smaller or larger than the lower limit of the confidence interval for the 
distance estimate, dmp = d-ri. 

The cell back radius Rb 'is not strictly needed, but provides enough information to 
improve the reliability of the method. Therefore in further embodiments of the 
present Invention the value of the back radius Rb Is set to zero. 

If diNF is larger than the serving cell back radius, Rb, the MS Is sufficiently far 
away from the serving BTS and is probably located outside of the back radius 
region. In this case, only the portion of cell in the direction of main radiation lobe 
of the serving BTS antenna is included in the confidence region. 

The inner radius of the confidence region is set to Ri = d-ri. However Ri = d-ri = 
dTA - X(Hyy2can be negative (i.e., if y«1 and the TA error statistics have long tails 
X(HYy2 can be larger than dxA). and If Ri is calculated to be negative, Ri should be 
set to zero. Moreover, when TA=0 It Is sensible to define the confidence region as 
a circle centred on the serving BTS coordinates; In which case Ri is set to zero. 
In summary, Ri Is defined as follows: 




max 



if TA = 0 
{04 -n) ifTA>0 



(197) 
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10 



The uncertainty radius of the confidence region (R2) is set equal to the width of the 
confidence interval of the distance estimate, n + rs. However, as a consequence of 
the adjustments made in (197) to determine Ri. the correct definition for R2 is 

R2 = d + rs-Ri. (198) 

The orientation angle of the confidence region (a) is determined by the orientation 
of the serving sector and the sector angular width. The width of the confidence 
region p is twice the serving sector width: 



{ 



a = ^g(rf) - A<f>s 

^ = 2A^5, (199) 



where A^^ Is defined to take into account conservatively the eventual dependence 
of the cell width on the distance from the serving BTS, d: 



15 The confidence region lias an angular width equal to twice the angular width of 
the serving cell. The sectorization of the ceil is therefore taken into account. The 
underlying reason for this is the assumption that all the mobile devices in 
communication with the BTS of interest are located in an arc defined angularly 
-A(t), <\|/(x,y)-<I),(d)< A(|>, (with an exception covered by equation (189) when Rb = 

20 0). This also produces the result that the confidence coefficient used to determine 
the circular crown, y, is equal to the confidence coefficient ^ used to define 
distances Ri and R2 



If diNF is smaller than the serving cell back radius, Rb, there is a large probability 
25 that the actual location of the MS lies in the back radius region. In this case the 
location estimate is calculated as previously but the confidence region is defined 
In such a way to include the whole back radius region. This means that the width 
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of the confidence region p is 2% and the orientation angle a provides redundant 
information as the confidence arc is defined as the whole circle (a is set to zero). 
Furthermore, the inner and the uncertainty radii are given as follows: 



In practical terms, the definitions above simply state that when the MS is very 
close to the serving BTS site, the effects of the back radius region are considered 
by treating the serving cell as if it was omni-directional. 

If the serving cell is omni-directional the concepts of sector orientation and 
angular width are meaningless. The best location estimate Is provided by using 
the co-ordinates given by the serving BTS site: 



where the confidence region is defined by a circle with radius equal to the upper 
limit of the confidence interval of the distance estimate: 




(201) 




(202) 




(203) 



d and rs are calculated as in the case of sectorized cell using equations (192) 
and (195). 



The calculations described above make use of statistical information from the TA 
measurement error. If the TA measurement error statistical information is not 
available, the definitions for location estimate and confidence region given above 
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Still apply; however the confidence coefficient is meaningless and d, n and rs can 
be defined only on the basis of the TA quantization rule: 

jx (cTb /2) s 138m if TA = 0 (204) 



(cTb/2)=277m if TA>0 



5 The algorithm presented above estimates the IVIS coordinates and the 
confidence region when both CI and TA are available. If no TA information Is 
available, the MS location estimate and its confidence region can still be 
determined by using the information carried by the CI only. In particular, the 
confidence region can be defined if an estimate on the cell radius, Rf is available 

10 while the location estimate can be provided even if Rf is not known. 

Where there Is no TA Information available, embodiments of the present 
invention examine the serving BTS to determine whether the serving cell is 
sectorized or not: 

15 

If the cell is omnidirectional and the TA is not available the location estimate Is 
according to embodiments of the Invention taken to be at the BTS coordinates: 

f X = Xs 

\y = ys (205) 

20 The confidence region of this estimate Is a circle centred on the BTS site. The 

radius of such region is obtained by scaling Rf by a factor -J^ . The factor Is 

chosen in order that a fraction B, of the total area of the circular cell of radius Rf 
is included In the confidence region: 
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10 



Xo = Xs 

a = Any value 



Ra = V?i?F (206) 



If the cell is sectorized and the TA Is not available, two alternatives exist for the 
location estimate. 

The first alternate as used in embodiments of the invention select the MS 
location estimate to be at the serving BTS site. In this estimate no cell radius 
(Rf) Information Is needed. 



X = Xs 

= ys (207) 



If the Rf is given, embodiments of the present invention calculate the estimated 
MS location as the mass centre of a simplified cell with the same shape as 
represented in figure 10 but with the back radius and the front radius scaled by a 
factor . to assure that only a fraction ^ of the total area of the original cell is 
15 considered. 

The coordinates of the serving cell's mass centre are defined as follows: 

~ L^d^dy ; Vmc = . / ydxdy 

M{S)Js M{S)Js ^208) 

where S is the border of the home cell, as explained in equation (189), and M(S) 
20 its area. By neglecting the dependence of the sector orientation and width on the 



distance and assuming a constant sector's angular width A^^ defined as 

A^ = ^m^^A<t>,(d) (209) 
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It can be shown that M(S) = Rf^A^ - (tc - A(j)7)RB^ and the Integrals in (208) 

solved In polar coordinates give the following expressions for the location 
estimate x = x,^ and y = y^ 

* 

X 

y 

V 

5 

The confidence region has the same shape as shown in figure 8, with the Inner 
radius equal to zero (R^ = 0). The origin of the confidence region (Xo,yo) is not at 
the BTS coordinates (Xs,ys) but is shifted along the axis defined by the sector 
orientation (in the direction opposite to the front side of the cell) by a distance 

10 Re'from the BTS coordinates. By using this definition of the confidence region, the 
back radius region is included (at least partially). The angular width of the 
confidence region ,p, is defined as twice the angle A^s^A^s. which is calculated 
In light of the different location of confidence region's origin and serving BTS 
coordinates. The confidence region's orientation, a, is defined according to the 

15 cell orientation, <|)s,'and the new variable A^s'- 

^0 = xs — Rb' coB<ps ; yo = ys — Rb' sin (f>s 
a = (f>s - A<l>s' ; fi = 2A(l>s' 

Ri = 0 ] Ba^ViiRF + RB') (211) 

Rb' and A^s are determined from the geometry as shown in Illustrated In figure 16. 

20 Figure 16 shows the geometry used in embodiments of the present invention for 
calculating the confidence region. The figure comprises a circle 1603, a first 
triangle 1601, a second triangle 1605, and a circular segment 1607. The circle has 
an origin, 1611, at the BTS location, Xs,ys, and has a radius Rb. The first 



+ 3 H>A?^-(.:S^)R,» COS <f>s 
, 2 (Rf^-Rb^) sin . , 



(210) 
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comprises a first vertex 1609 at location Xo.yo, a first side 1623, connected at one 
end to tlie first vertex, 1609 wtiicli passes the origin of tlie circle, Xs.ys, a second 
side 1621, connected at one end to the first vertex and arranged so that the angle 
at the first vertex 1609 defined by the two sides 1623,1621 is A^s- The vertex 
5 1609 of the first triangle is located a distance Rb' from the origin of the circle. The 
second triangle 1605 lies within the area of the first triangle 1601, and comprises 
a first vertex located at the origin of the circle, a first side 1629 connected at one 
end to the second triangle's first vertex and forming part of the first triangle's first 
side, and a second side 1627 connected at one end to the first vertex and 

10 arranged to have an angle between the first side 1629 and the second side 1627 
of A<|), . Both the first triangle 1601 and the second triangle 1605. also comprise a 
common right angle vertex 1631, and a common side 1625. The common side 
being arranged as the line perpendicular to both triangles' first sides and 
connected to the point where the first and second triangle's sides intersect 1635. 

15 The segment 1607 comprises two radii of length Rf with an origin at the origin of 
the circle 1603 , Xg.ys, and an arc defined between the ends of the radii. The first 
radius is common to the second triangle's second side 1605, and the second 
radius is an angle 2 A(j)j from the first. Using simple trigonometry, the height of the 
first and second triangle's second side can therefore be shown to be equal to 

20 RFSln(A<|), ). Figure 16 further has a line 1651 connected to and perpendicular to 
the first triangle's first side 1623 at the circle's origin 1611, and a second end at 
the intersection of the line and the first triangle's second side. The line 1651 has a 
length h. 

25 Depending on whether is larger or smaller than n/2, embodiments of the 
present invention determine Rb' and A^s' as explained below. 
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10 



15 



20 



If Rf » Rb and A^^ > n/2, it can be seen that only a very small area of the circle 
defined by Rb is left out of the confidence region if Re' is defined as being equal to 



Rr- 



Rb' = 



(212) 



Using this definition for Rb' the analytical expression for A^s' is calculated from the 
geometry shown in figure 16. The result is: 



^4>s' = TT — cos ^ 



■Rf cos(7r - A^^) - Rg' 



VCiJjr cos(7r - S^) - Rb'^ - Rf^ sm2(7r - S^. ' 



(213) 



Using the geometry shown In figure 16 A^s' is obtained as follows: 

A<l>s' = tan~^ 



■Rj-sin(A^g) - h 



Rf cos(A0s) 



and using this result Rb' is easy to calculate as follows: 

h 



(214) 



Rb' 



tan(A^.9') (215) 



The length h can be approximated as h « Rb. This approximation is only performed 
after checking that the confidence region is too small and does not Include a large 
enough part of the circle with radius of Re (i.e. the value of Rb' Is too small) or 
conversely that Rb' is too large leading to a too large confidence region. To check 
for these cases two chosen parameters are introduced: Smin (taking care of the too 
small Rb' example - where 1.5 is a typical choice for 5min) and 5max (taking care of 
the too large Rb' example - there 4-5 is a typical choice for 6max). 
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If Rb' turns out to be too small, I.e. If 8minRB > Re', then Rb' Is re-deflned In 
embodiments of the Invention as 

Rb' = <5min-RB 



(216) 



and A<l)s' is obtained using the geometry as shown in Figure 16: 



Rp sin(A^5) 



Rp cos{A(f>s) + Rb'I (217) 



If Rb' turns out to be too large, i.e. if SmaxRe > Rb' then Rb' is re-defined as 

Rb^ ^ SjobxRe (218) 



and s' is calculated from equation (217). 



The proposed location procedure incorporating at least some of the location 
methods utilizing CI,TA or andRXinformation is shown in figure 17. The procedure 
can incorporate location methods presented above (CI, CI+TA and CI+RX 
algorithms). The process of which is described as follows: Firstly, in step SI the 
timing advance and received signal levels measured by the mobile station from the 
serving and neighbour cells are collected. 



Next, the relevant algorithm input parameters and radio network parameters such as 
the base transceiver station coordinates for the measured cells are collected in step 
S2. 



In step S3, the measurements and the network data are analysed to select which 
measurements have to be used in the location calculation. Cells that do not 
contribute significantly to improving the location accuracy are removed from the 
measurement set. As a result of the cell selection procedure, serving cell information 
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can potentially be removed and/or neighbour cell infomiatlon can possibly be totally 
or partially removed. 

In step S4, a cell identity/cell identity plus timing advance location estimate is made. 
5 If serving cell information is available, cell identity or cell identity and timing advance 
based location algorithms such as described hereinbefore can be used to determine 
a location estimate and/or a confidence region for it. The resulting location estimate 
is called in this document, the "CI/CI + TA Location Estimate". 

10 In step S5, an RX location estimate is made. If the neighbour cell Information is 
available, a CI + RX based location algorithm is applied. This may be any one of 
algorithms A to F described hereinbefore or indeed any other suitable algorithm. This 
is used to estimate the mobile station location. In selecting the CI + RX based 
algorithms, the following criteria may be taken into account: 
15 • cell sectorisation. For omnidirectional cells, CI + RX based algorithms can be 
used. For example, algorithms D (equation 154), E (equation 156) or F (equation 
158) may be used. 

For sector cell cells, CI + RX based algorithms such as algorithms A, B and C 
20 described hereinbefore can be used directly. Alternatively, algorithms D, E and F can 
be extended to the case of sector cells as described previously. 

• Absolute/relative level observations 

25 Algorithm C described hereinbefore makes use of relative level observations, that is 
the difference in received signal level measurements from pairs of base transceiver 
stations. It should be appreciated that algorithms A, B, D, E and F use absolute level 
observations. 
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• Closed form/iterative algorithms 

Algorithms A, B, C, D and E are iterative algorithms. Algorithm F is a closed form 
algorithm. It determines the mobile station coordinates as a weighted average of the 
5 coordinates of the base transceiver stations from which the received signal level 
measurements have been collected. One analytical expression (equations 163 and 
164) and three empirical approximations for the weights are also provided as 
discussed above 

• RX-based domain determination 

10 

The RX-based algorithms may require determination of the "domain" in which the 
RX-based location estimated is constrained to lie. Such a domain can be defined 
using, for example, any of the definitions provided hereinbefore. 

15 The obtained location estimate is called for convenience "RX location estimate". 

In step SB an RX direction location estimate is obtained. If serving cell information 
and the "RX location estimate" are available, the direction from the serving base 
transceiver station to the "RX location estimate" is calculated. This direction is used, 
20 instead of the actual antenna direction, as an input for a CI + TA based algorithm, as 
explained above. The resulting location estimate is called, for convenience, "RX 
direction location estimate". 

In step S7, a virtual BTS location estimate is calculated. 

25 

If serving cell information and "RX-direction location estimate" are available, the 
coordinates of the "RX-direction location estimate" is used as an additional (omni- 
directional) neighbour cell for the RX based algorithm. This additional neighbour cell 
and the RX level measurement associated with it can be called the "virtual base 
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transceiver station" and "virtual RX level" respectively. The value for the virtual RX 
level measurement can be selected for example as the maximum RX level from the 
measured neighbours or as the average RX level from the measured neighbours. 
The resulting location estimate is called "the virtual BTS location estimate". 

It should be appreciated that a location estimate obtained by any of the location 
methods discussed can be used as an additional or virtual base station 
measurement. The virtual base station measurement is associated with a virtual 
measurement, for example the virtual RX level, and the set of real measurements 
and the virtual measurement is reprocessed using any of the location methods 
described. 

In step S8, any one of the obtained location estimates: CI/CI + TA location estimate; 
RX location estimate; RX-direction location estimate; RX location estimate; RX- 
15 direction location estimate; or virtual BTS location estimate can be selected as the 
location estimate delivered by the location process. 

Embodiments of the present invention are particularly flexible. For example, if the TA 
measurements and/or serving cell information are not available, the process is still 
20 able to provide location Infomnatlon - such as the RX location estimate. Likewise, if 
neighbour information and/or RX level measurements are not available, the process 
is still capable of delivering location estimates using the CI/CI + TA location estimate. 

Embodiments of the present invention enable both the RX level and TA 
25 measurements to be used effectively. This Infomnatlon is currently provided in the 
standard GSIVI networks. The combined use of timing advance and RX level 
measurements "the virtual BTS location estimate" gives the most accurate location 
estimates in the regions of the highest density of the mobile users. In embodiments 
of the present invention, one of the estimates is selected. Where more than one 
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estimate can be obtained, an averaging or alternatively a weigfiting of the results may 
be provided. 

Wliere a number of different location estimates can be determined, only one of the 
methods may be used. Any suitable criteria can be used for determining which of the 
methods is to be used. 

Whilst embodiments of the present invention have been described particularly in the 
context of an arrangement where a number of different location methods are used, it 
should be appreciated that in some embodiments different ones of the location 
algorithms described hereinbefore can be used alone. 
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CLAIMS 

1 . A method of estimating the location of a mobile device, comprising the steps 
of: 

5 collecting location infomiation; 

selecting at least one of a plurality of different location methods to provide a 
location estimate; and 

providing a location estimate based on the at least one selected location 
method. 

10 

2. A method as claimed in claim 1 wherein said at least one location method 
comprises at least one of the following methods: 

a method using cell identity information; 

a method using cell identity information and received signal strength; 
15 a method using cell identity information and timing advance information; and 

a method using cell identity information, received signal strength information 
and timing advance information. 

3. A method as claimed in claim 1 or 2, comprising the step of determining a 
20 virtual base station estimate. 

4. A method as claimed in claim 3 when appended to claim 2, wherein said 
virtual base station estimate is determined using at least one of the methods of claim 
2. 

25 

5. A method as claimed in claim 3 or 4, wherein said virtual base station location 
estimate coupled with at least one virtual measurement and at least one real 
measurement and said at least one virtual measurement is processed using a 
location method. 
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6. A method as claimed in claim 5, wherein the at least one real and the at least 
one virtual measurements are processed using a location method as defined in claim 
2. 

5 

7. A method as claimed in claim 5 or 6, wherein a value for the virtual 
measurement is one of measured levels, a combination of measured levels, and an 
average of measured levels. 

10 8. A method as claimed In any preceding claim, wherein said at least one 
location method is selected in dependence on the location information available. 

9. A method as claimed in any preceding claim, wherein a plurality of location 
estimates are determined and at least one estimate is used to provide said location 

15 estimate. 

10. A method as claimed in any preceding claim, wherein said location information 
is collected by said mobile device. 

20 11. A method as claimed in claim 10, wherein said mobile device is arranged to 
measure a level of at least one type of information. 

12. A method as claimed In any preceding claim, wherein said location information 
comprises at least one of timing advance information and received signal level. 

25 

13. A method as claimed in claim 12, wherein said received signal level is an 
absolute received signal level or relative received signal level. 
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14, A method as claimed in any preceding claim, wherein said mobile device Is in 
a cellular communications device. 



15. A method as claimed In claim 14, wherein said information is collected for a 
5 serving cell of the mobile device. 

16. A method as claimed in claim 14 or 154, wherein said information Is collected 
for at least one neighbouring cell. 

10 17. A method as claimed in any of claims 14 to 16, comprising the step of 
selecting the or each cell in respect of which location Information is collected. 



18. A method as claimed in any preceding claim, wherein a location estimate is 
provided using the following algorithm 
15 Calculate the total attenuation experienced by a signal transmitted by the 

i-th BTS while propagating toward a mobile station where i-th level observation 
is L') by subtracting from the i-th measured received power, P\, the maximum 
power radiated by the i-th BTS. PVmax: 



Stack the level observations from N BTS's in vector L: 

(12) 



25 



Solve the minimization problem: 



^ 2 

X 



= axg 



min F {x,y\ aj^^ 



X 



(13) 
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Where the cost function F (x,y; aa^) is defined as follows: 

F {x,yi <7„2) = ln<r„2 + La \Ti.ix,y)\ + [L - mi,(a;,y)]^ri,-i(x,y) [L - mi.(a;,y)] 

(14) 

5 and 

mi,(a;,y) = [fiiM, . . . ,/iif(a;,y)]^ ^^^^ 
lAM = -PL* {d*{x,y)) - API, (iP'ix,y)) ^^^^ 



10 



(17) 



19. A method as claimed in any preceding claim, wherein a location estimate is 
provided using the following algorithm 

Calculate the total attenuation experienced by a signal transmitted by the 
l-th BTS while propagating toward a mobile station where the i-th level 
15 observation is L' by subtracting from the i-th measured received power, Pr, the 
maximum power radiated by the i-th BTS, P^max: 



(18) 



Stacl< level observations from N BTS's in vector L: 
20 L= ^^gj 



Solve the minimization problem: 
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X 

y 



— arg ^ nnn F 
X ^ 



y 



(20) 



where the cost function F(x,y) is defined as follows: 



(21) 



and Dxyis the domain of existence of x and y. 
Calculate al as 



(f„2 = F {x,y) (22) 



10 20. A method as claimed in any preceding claim, wherein a location estimate is 
provided using the following algorithm: 

Calculate the total attenuation experienced by a signal transmitted by the i- 
th BTS while propagating toward a mobile station where the i-th level 
observation Is L') by subtracting from the i-th measured received power, Pt, the 

15 maximum power radiated by the I-th BTS, P't, 



max- 



(23) 



20 



Calculate the j-th level difference observation by subtracting the j-th level 
observation from the level observation taken as reference: 



Stack the N - 1 difference of level observations in a vector D: 

(25) 



D=[p^...,i?^]' 
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Solve the minimization problem 



X 



= arg ^ min F {x,y) 

X 



y 



(26) 



where 



(27) 



and 



/4>(^.y) = - [PL^ (dHa;,y)) - Pi/ (ti'Caj.y))] - [^P^ {^\x,y)) - API {f{x,y))] 



(28) 



10 Dxy is the domain of existence of x and y. 

21. A metliod as claimed in any preceding claim, wiierein a location estimate is 
provided using an algorithm solving the following equation in x and y: 



15 



t=l 

f:^W)(2/-!/0=O 



where 



{B^ logio (x,y) - z< + ^ 



2£7*=^llllO 



] 
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22. A method as claimed in any preceding claim, wherein a location estimate is 
provided using an algorithm solving the following equation in x and y: 



N 

E 

ts=l 



N r 

L 1=1 



|R| 



; {x,y)&v 



23. A method as daimed in any preceding claim, wherein a location estimate is 
provided using an algorithm based on the following equation: 



X 



24. A method as claimed in any preceding claim, wherein said location estimate is 
10 provided by one of a iterative and a closed form method. 

25. A method as claimed in any preceding claim, wherein said location estimate is 
provided by one of a linear and non linear method. 

15 26. A system for estimating the location of a mobile device, comprising: 
means for collecting location information; 

means for selecting at least one of a plurality of different location methods to 
provide a location estimate; and 

means for providing a location estimate based on the at least one selected 
20 location method. 
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Fig. 1 
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Fig. 4 

Example of antenna gain 
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Fig. 5 
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Fig. 8 

D from Cell Identity (and Timing Advance) of serving cell 
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Fig. 10 

Anaiytical definition for Dj 



1007 




SUBSTITUTE SHEET (RULE 26) 




SUBSTITUTE SHEET (RULE 26) 



wo 2004/023155 




PCT/ra2002/003644 



9/12 

Fig. 12 Plots oflio in equation (45) as a 
function of patti-loss/attenuation 

IfO as a function ofz' ai = 8, A'= 125, B'=40 
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Fig. 14 

The constant'TA circular crown defined in equation (54) 
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Fig. 15 

Definition of estimate of distance and its confidence interval 
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Fig. 17 
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